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ABSTRACT 


German Army light helicopter transportation regiments operate 
45 Bell UH-1D helicopters to support demanding missions throughout 
Europe. Maintenance period scheduling, major exercise and regular 
mission assignment decisions directly influence the readiness of 
the helicopter fleet. Currently, all planning is done manually, 
which is unstructured and time consuming. This thesis describes a 
decision support system designed to assist with maintenance 
planning and mission assignment. The yearly maintenance and event 
scheduling problem and the short term mission assignment tasks are 
formulated and solved as elastic mixed integer linear programs. 
Resulting yearly schedules and short term sortie plans are both 
generated in a fraction of the time previously required with 


solution quality superior to their manual counterparts. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in 
this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
available, to ensure that the programs are free of computational 
and logic errors, they cannot be considered validated. Any 
application of these programs without additional verification is at 


the risk of the user. 
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EXECUTIVE SUMMARY 


This thesis develops two optimization models and proposes 
an interface for a self sufficient personal computer based 
decision support system as an interactive instrument to 
construct reliable and completely organized helicopter usaqe 
and maintenance plans. The optimization models assist with 
yearly maintenance and event scheduling, and short term 
helicopter - mission assignments. Computational experience 
shows both yearly schedules and short term sortie plans are 
generated in a fraction of the time previously required with 
solution quality superior to their manual counterparts. 

These results are for a German Army light helicopter 
transportation regiment operating 45 Bell UH-1D helicopters. 
For such regiments, maintenance period scheduling, major 
exercise, and regular mission assignment decisions directly 
influence the readiness of the helicopter fleet. The planning 
supervisor in the regiment's maintenance and repair battalion 
strives to keep high technical and operational standards while 
meeting all necessary inspections (16 different inspections 
during a 1200 flight hours cycle), satisfying all mission and 
exercise requirements, equitably using the helicopters, and 
smoothly operating the maintenance facilities. Done manually, 
these tasks are unstructured and time consuming. 

The two optimization models provide valuable assistance to 


the planning supervisor. The yearly planning model assigns 


УЛТ 


helicopters to inspections and to exercises while observing 
monthly planned flight hours and operational guidelines (а 
desired level of flight hour reserve, a percentage range of 
operationally ready helicopters, an upper level of monthly 
flight hours per helicopter, and inspection capacities). The 
short term model assigns helicopters to missions while 
observing the technical status of each helicopter (remaining 
flight hours to next inspection, planned time of next 
inspection, availability) md ul filling ειν mission 
requirements with respect to flight hours, and equipment. 

Resulting yearly schedules and computed short term mission 
assignment plans are face-valid (i.e. judged implementable by 
expert opinion), superior to their manual counterparts (i.e. 
fewer planning conflicts) and generated in a small fraction of 
the time previously required. 

Included in this study are the system's structure, the 
development and implementation of the two optimization models, 
computational results, the interface proposal, and 


explanations regarding the planning process and intended user. 


J 





г ІШІКОрУСТТОМ 


German Army light transport helicopter regiments operate 
45 Bell UH-1D helicopters (see Figure 1) in support of a corps 


with three army divisions. A maintenance officer or "Leiter 





Figure 1 


Einsatz" (the author's assignment from 1989 to 1991) in the 
regiment's maintenance and repair battalion supervises the 
helicopter fleet maintenance planning and mission assignment. 
He strives to keep high technical and operational standards 
while: 

ӨФӨ Meeting all necessary inspections, 

® Satisfying all mission requirements, 

ө Fulfilling special events such as NATO exercises, 

® Equitably using the helicopters, 

ө Smoothly operating all maintenance facilities. 
This thesis derives and solves integer linear programs to 


assist with maintenance planning and mission assignment. 


A. BELL UH-1D INSPECTION SYSTEM 

The Bell UH-1D maintenance cycle (in its German version) 
consists of 1200 flight hours and contains 16 а реа 
inspections in two levels (C and D). The C-level (see Table 
1) includes relatively easy maintenance and part replacement, 
the D-level (see Table 2) contains costly repair overhauls 
lasting three to five weeks. Each inspection takes place 75 


flight hours after its predecessor in the following order: 


C1,C2,C1,D1,C1, C2, C1 ο Ci, ieee СТ Co CT I E 


C and D inspections are accomplished in disjoint facilities 
with different technical personnel. The C inspections, being 
relatively easy, can be performed without using fixed 


installations (e.g. during exercises). 


TABLE 1 







TABLE 2 


D-LEVEL (Repair) 
DURATION Ε.Κ 


; 


B. CURRENT PLANNING PROCESS 










The maintenance officer in charge of the planning and 
mission assignment tasks has a staff of up to six soldiers 
experienced in the fields of aircraft repair and maintenance. 
Despite access to a personal computer (PC), maintenance and 
mission assignments are currently based on poorly documented 
manual procedures which rely on complicated charts and 
overview boards. 

The planning process breaks into two related pieces, 
yearly and short term. The yearly maintenance and event 
schedule takes a three man team up to five weeks to produce. 
The final product contains the flight hours assigned to each 
helicopter each month, helicopters assigned C and D 


inspections each month, and helicopters assigned to fly 


Special events (e.g. operations of the Allied Mobile Force in 
northern Norway or in eastern Turkey, which require up to 12 
helicopters to fly a"totat" "orf morethan ОООО The 
short term plan provides individual helicopter mission 
assignments for up to one week. Time needed to develop the 
short term plan varies substantially depending on mission 
requirements and helicopter availability. A typical short 
term plan requires one man approximately two hours. 
l. Yearly Maintenance and Event Schedule 
Input to the yearly maintenance and event schedule 

includes: 

e The total number of hours the helicopter fleet should fly, 

€ Special event requirements, 


e Pilot Combat Training Programs (CTP, provide required 
monthly instrWet ionad erage sours), 


€ Operating data from past years on availability апа 
reliability of the (aging) helicopter. 


The yearly schedule has four primary operational 
guidelines: 
e The flight hour reserve, 
ө Тһе percentage of operationally ready helicopters, 
ө Тһе upper level of monthly flight hours per helicopter, 
ө Тһе maximum monthly C and D inspection capacities. 
The flight hour reserve is each helicopter's available flight 


hours to the next D inspection, summed over all helicopters. 


ІШ ісеса! level of 6,750; hours (number of helicopters * 0.5 * 
hours between D inspections or 45 * 0.5 * 300 = 6,750) has 
historically been a good planning factor. A level above 7,500 
hours has endangered future equitable use of the maintenance 
facilities (i.e. a disproportional number of helicopters 
require imminent inspections), a level below 5,500 hours has 
endangered the capability of the regiment to fulfill all 
required missions. 

The percentage of operationally ready helicopters measures 
the number of helicopters not in a D inspection and with 
remaining hours to the next D inspections greater than zero. 
A level between 70% and 90% is desired. 

No more than 30 flight hours should be assigned to each 
helicopter each month, but violations are sometimes necessary 
(e.g. for events). 

mie megmal output capacity for D inspections is three per 
month. If planned well ahead, an output of four per month is 
achievable, but reserve capacity for exception repairs is lost 
and an equally high output in the succeeding month is 
unlikely. The C-inspection level capacity is easier to 
manage. Output variations from one to six helicopters in one 


week are possible. 


2. Short Term Helicopter - Mission Assignments 
The short term planning process attempts to pick the 
correct helicopter for each mission while observing: 


e The number of flight hours each helicopter has until the 
next inspection 


e The planned time and importance of the next inspection, 
Φ The mission's flight hour requirement, 


e The current equipment status and the mission's equipment 
requirement, 


e The possible multiple use of the helicopter for non- 
simultaneous missions. 


This daily task requires experience and talent. 
Unfortunately, an organized method of meeting the listed 


criteria is often blocked by time considerations. 


C. OBJECTIVES OF THIS THESIS 

The objective is to develop underlining algorithms and 
propose an interface for a self sufficient PC based decision 
support system as an interactive instrument to reach reliable 
and completely organized helicopter usage and maintenance 
plans. 

The approach undertaken in this thesis is practical and 
centered around two optimization models. 

Chapter II addresses similar concepts in the existing 
literature. Chapter III contains the developed algorithms, 
the proposed interface and both user and planning session 


descriptions for the decision support system. Chapter IV 







is the computational performance of the developed models 


.er V provides conclusions. 


II. RELATED RESEARCH 


A literature search did not reveal a model with the 
capabilities needed by the German Army Aviation maintenance 
officer. His tasks are very specific, dependent on the German 
version of the Bell UH-1D preventive maintenance system and 
the local maintenance resources. The combination of a yearly 
plan with a short term mission assignment system, the 
necessity of using this system independently and the exclusion 
of manpower, budget and logistic issues make the proposed 
model unique. 

However, decision support and expert systems for aviation 
maintenance activities have been studied for other aircraft 
and organizations. Hackett and Pennartz (1982) provide the 
basis of a decision support system for the maintenance 
aircraft scheduling process of an United States (US) Air Force 
B-52 wing. They start with the principle that computerization 
without proper management does not improve overall 
performance. They therefore first establish maintenance 
responsibilities and management procedures for an operational 
cycle that considers everything from the yearly flight program 
to the daily scheduling update operations. They characterize 
maintenance scheduling as a complex process with a high degree 


of uncertainty, strict requirements, binding constraints and 


mew ficient. guidance. In their view, successful usage of 
Optimization models in the different planning stages depends 
on appropriately reducing the complexity of the problem with 
simplifying and structuring techniques. They do not, however, 
develop optimization models or algorithms. Key aspects of a 
decision support system are reported as: 

® Communicability with the user, 

èe Robustness even for extreme cases, 


e Ease of Control. 


Shenolikar (1983) describes a decision support system for 
automatic test equipment systems operations management, which 
is closely related to aircraft maintenance affairs. 
Optimization models are again addressed, but not formulated. 
The key elements of a generic decision support system are 
developed as follows: 


e Knowledge base (with proposed algorithms and solution 
models), 


e Data directory and data base, 
u porc generator, 


@ Communications (interface) manaqer. 


A series of theses (see Christensen and Pasadilla,1991) 
advised by Professor Martin J. McCaffrey (Department of 
Administrative Science, United States Naval Postgraduate 


School) develops a Naval Aviation Maintenance Organizational 


Activity Strategic Information System (OASIS) and an Expert 
System Advisor for Aircraft Maintenance Scheduling (ESAAMS). 
These complex systems address flight and maintenance 
activities, but include also human resources (manpower 
management, training and qualification), monetary management 
(budget considerations and accounting) and inventory 
structures for logistical support. These programs deal with 
the combined use of present information systems of the US 
Naval Air Systems Command and design applicable knowledge 
bases, databases and graphical interfaces. Again, these 
systems only propose optimization models and algorithm as 
future work. 

In addition to these expert systems for planning aircraft 
related maintenance and use, the literature search did reveal 
an optimization model concerned with helicopters and one 
dealing with military aircraft sortie planning. The "Phoenix" 
model (Brown, Clemence, Teufert and Wood, 1991) includes 
procurement and retirement schedules for the US Army’s 
helicopter fleet, handling 16 different helicopter types over 
25 years. The mission assignment problem is described by 
Wallace (1992), who develops sortie optimization tools for the 
US Air Force. 

The search for related maintenance planning systems in use 
by commercial airlines was not successful. Commercial 
airlines are concerned with a different set of constraints 


(costs) and resources (see Talluri and Gopalan, 1993). 
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III. CONCEPT OF THE DECISION SUPPORT SYSTEM 


A. STRUCTURE AND OVERVIEW 
The system organization (see Figure 2) resembles the 
underlying manual planning process and consists of two almost 
independent parts: 
e The yearly maintenance and event planning system, 


e The short term planning system for mission assignment 
covering up to one week. 


Planning Data 


D 


Yearly Planning 
System 


Update Process 
User Preferences 
Mission Data 


5 Schedule „Ф д 


Maintenance 


| | 
Schedule | Short Term 
Event Decision 


Flight Hour Planning System 
Distribution 
Statistics | 
Graphical p 
Overview | | 
| Mission Assignments | 


— 





Figure 2 
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Each subsystem has the ability to respond to user preferences 
and demands. The scheduled inspection periods are input to 
the short term planning system which forms a connection 
between the two parts. There is no direct feedback out of the 
short term planning process into the yearly system, because a 
continuous change of the yearly plan destroys the underlying 
planning policy. There are good possibilities to get a yearly 
plan "back on track" when unforeseen trouble strikes. The 
experienced user has to decide, when a complete renewal of the 
yearly plan is necessary. Every subsystem result can and 
should be adapted manually. Results are proposals, providing 
feasible starting points for further planning and adapting to 
reality. All the important existing constraints are present 
in the different models, but elastic violations are possible 


(like real decision making) by paying adequate penalties. 


B. THE YEARLY PLANNING SYSTEM 
l. Basic Specifications 

The yearly planning system forms the basis for all 
maintenance scheduling and mission assignment. It recognizes 
initial conditions and requirements from the planning data and 
implements the planning policy. The essence of this subsystem 
is an integer linear program. The necessary input can be 
expressed as: 


e Planned flight hours per month for the fleet, 


82 


Special event requirements (time, flight hours, number of 
helicopters, maintenance possibilities), 


D-level maintenance capacity, 


Status of each helicopter (remaining flight hours to next 
D inspection). 


final output includes 
Recommended flight hours for each helicopter each month, 


Recommended D inspection decision for each helicopter each 
month, 


Helicopter assignment to special events, 


Monthly statistics with respect to flight hour supply and 
helicopter availability. 


A Pascal program transforms the yearly solution (containing 


monthly decisions) into a solution containing weekly results 


and C inspection decisions. A description of this program 


follows the yearly planning model. 


2. The Yearly Planning Model 


The model specifications in a basic format can be 


G@escribed as follows: 


Be INDICES: 
Г ШЕШЕН (е.а. 1,2,..912), 
h Ae oper identification (e.g 1,2,...,45), 
e event (e.g. 1,2,...,E). 
W DATA: 
Pinko, planned flight hours for month t, 
EVENTHRS, flight hours required for each helicopter 


participating in event e, 


142 


EVENTNMB, 


DMAX 


PRODHRS 


ОРАР 


MININSP, 


MAXHRS 


B VARIABLES: 


number of helicopters required for event 
e, 

maximum number of D inspections per month, 
additional flight hon (e.g. 300) 
obtained per inspection, 

desired flight hour supply (e.g. 6,750), 
minimum number of inspections in month t, 
maximum regular monthly flight hours for 


each helicopter (e.g. 30). 


BINARY VARIABLES 


Zh 


"e 


ЕЕ 


one if helicopter h is assigned to event e, 
one if helicopter h is assigned a D inspection 


in month t 


POSITIVE VARIABLES 


Χμε 


A P 


БЕМНЕ5,. 


B MODEL: 
* Minimize the 


hours. 


flight hours assigned to helicopter h below 
MAXHRS in month t, 

flight hours assigned to helicopter h in month 
t exceeding MAXHRS or in addition to an 
inspection or event assignment, 

flight hours until next D inspection for 


helicopter h at end of month t, 


sum of all weighted elastic penalties and XE,, 
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w bece to the-constraints: 


ШЕ) 


ο) 


45 Ε 

РОО XE t EVENTHRS,*Z, .] *PLHRS, curte 
h=1 e=1 

45 

Y, Y, < *DMAX Sele 
h=1 

45 

De ο ΜΜ, E=1..12 
h=1 


ο απ -REMHRS, .,*tPRODHRS*Y, .-X, .-XE, .- 
E 
y. EVENTHRS,*Z, , 


езі 


οσο ο 12 


45 


Y. КЕМНК5, ,><ОРТЗИУР t=1..12 
h=1 

12 E 12 

1 ΕΝΤΗΕΘ,:2, .|=()) PLHRS,) /45 

С-1 егі 


t=1 


5 


15 


45 
(73) Y. Z, ,"EVENTNMB, e-1..E 
h=1 


(8) 2, Үр, (X4, / MAXHRS) «1 hz1..45,e-1..55t-1 9 


m REMARKS: 


signify elastic constraints. 


e REMHRS,, provide the initial status of each helicopters: 


mE DESCRIPTION OF THE CONSTRAINTS: 

(1) Meet planned flight hours each month or incur an elastic 
penalty. 

(2) Comply with the maintenance capacity or incur an elastic 
penalty. 

(3) Perform a minimum number of inspections each month. 

(4) Calculate remaining hours until the next inspection for 
each helicopter at end of each month. 

(5) Provide the desired flight hour supply or incur an elastic 
penalty. 

(6) Provide equitable use of the helicopter fleet or incur an 
elastic penalty. 


(7) Enforce number of required helicopters for each event. 
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(8) Allow each helicopter to participate only in a special 
event, be inspected, or be assigned regular flight hours 
Χμε without penalties each month. Assignment of flight 
hours XE, e despite an event or inspection is possible, but 


incurs the penalty associated with the variable XE,,. 


С. TRANSFORMATION OF OPTIMIZATION RESULTS 

А PASCAL program named  "MOSCH" (Monthly SCHedule) 
transforms the model solution into an editable form (see 
Appendix C). The objective is to translate the monthly 
decision variables into weekly results. Instead of operating 
with a general maintenance duration time for the D inspections 
of one month, the more realistic values for the respective 
D1,D2,D1a and D3 overhauls are utilized. The three required 
C inspections between every two D inspections are also added 
(recall the relative ease of planning these inspections). The 
basic approach can be described as follows: 


€ Process each helicopter ordered by initial flight hours 
until next D inspection. 


e For each planned D-level inspection, schedule the 
inspection as soon as possible starting two weeks prior to 
the beginning of its planned inspection month. Ensure no 
planned hours or event conflict exists. 


e If the starting week is feasible, check inspection 
workload and choose next available week with acceptable 
workload level. 


e Schedule the inspection over its actual duration and check 
again for scheduling conflicts. 


e Schedule all intermediate C inspections based on remaining 
hours to the next inspection for each helicopter. 


T7 


Ө Note C inspection requirements during events. 
e Print the schedule, a graphical overview and all required 
statistics. 
D. THE SHORT TERM PLANNING SYSTEM 
1. Basic Specification 
The short term planning system assigns helicopters to 
missions while observing requirements of the yearly plan. It 
is proposed as an interactive program with a graphical user 
interface and an optimization model to perform mission 
assignment. Additional characterizations of the complete 
short term planning system are given later in this chapter. 
The following addresses the optimization model, which can 
be used independently from the proposed system. Necessary 
input requirements for this mission assignment model are: 


€ Status of each helicopter with respect to availability, 
equipment, flight hours and inspection plan, 


èe Mission requirements with respect to fiight hours and 
equipment. 


The output includes an assignment proposal for each mission. 
Some sorties require a spare or backup helicopter to ensure 
mission success. The optimization model also decides the 
assignment of spare helicopters. 
2. The Short Term Helicopter - Mission Assignment Model 
This model chooses the best helicopter and (if 
required) spare helicopter for each mission. It takes the 


following helicopter properties into account: 


Το 


A user preference for mission assignment, 


The relative importance of the next inspection, expressed 
in weeks required for the completion, 


The time gap until the next planned inspection, 

The flight hours remaining until the next inspection, 
The current equipment condition, 

Restrictions for night or instrument flights, 


The overall availability during the planning time frame 
(по ongoing mission, event assignment or inspection). 


Mission requirements are: 


The 


Beuxpment, 
ENESht hours, 
Backup demand, 
Time window, 
Et or instrument flight capability. 
equipment requirement is one of the following: 
ll seats, 
5 seats, 
400 kg exercise load, 
Internal tank left, 
Enternal tank right, 
Winch. 


number of items between any two equipment types on the 


above list (noncircular) represents the relative difficulty or 


time needed to change between the two equipment types. 


Therefore, a change between 400 kg exercise load and 5 seats 


is relative easy to do in comparison to a change between 11 


ШӘ 


seats and a winch. This approach is chosen, because data for 
a more realistic change-over matrix is currently not 
available. 

A mission - helicopter assignment is only feasible, if the 
helicopter is available (expressed in a 0-1 availability 
subset), the helicopter is not assigned to a simultaneous 
mission, and mission requirements don't collide with flight 
restrictions of the helicopter. The mission requirements fall 
into one of the following codes: 

@ 0 = no operation limitations, 

® 1 = helicopter restricted to Combat Training Flights (CTP) 
because of insufficient technical standard (e.g. non 
critical vibraciono) 

Ө 2 = helicopter restricted to Visual Flight Rules (VFR) 
missions because of absence of Instrumental Flight Rules 


(IFR) equipment, 


Ф З = helicopter restricted to daylight ҒІтаһра -ᾱ-- 
because of absence of night sight equipment. 


The model structure can be described as follows: 


RE INDICES: 
h helicopter identification (e.g. 1,2,...,45), 
m mission identification (eem 7 MT 
il simultaneous mission group (e.g. 1,2,...,I). 
E DATA: 
Reminsp, remaining hours to the next inspection for 
each helicopter h, 
Length, flight hours for mission m, 
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Costuse,, penalty for assigning helicopter h to mission 


m, 
Spare subset of all missions which require a spare 
helicopter, 
Overlap; set of all missions in group i. 
B VARIABLES: 
EM, z binary assignment decision for helicopter h 


and mission m, 
B a binary assignment decision for helicopter h 


used as a spare for mission m. 


B MODEL: 


арат: ле 


45 М 
ШИ Costuse, ,* (FM, ,+FS, ,) 
H-1 m-1 


EEubgject to the constraints 


(1) y FM, l m-1..M 
іі 
45 

(2) AES. ul m-1..M s.t. m € Spare 
іі 
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M 


(3) Y ` [Length,* (FM, ,*FS, ,) ] sReminsp, h=1..45 
m=1 

(4) З и ІЕМ, ,%Е5,,1<1 hw. .45989 τ 
M 

(5) 9” ЕМ, 1<51 h=1..45 
med 


mE DESCRIPTION ΟΕ ΤΗΕ ο ΙΝΕ: 


Тһе 


tos 


Assign exactly one helicopter to mission m. 

Assign a spare helicopter to mission m if required. 

Each helicopter's flight hour use can not exceed its 
remaining hours until inspection. 

A helicopter can not be assigned to more than one mission 
being conducted at the same time. 


Assign each helicopter to at most one mission or incur an 


elastic penalty. 


variable Costuse,, 1s a combination of weighted penalties 
not choosing a helicopter with: 


the highest priority score, 


the most expensive following inspection (expressed in 
required completion time), 


үй 


èe the least expensive equipment changeover (expressed in 
equipment code differences), 


e the smallest time gap to the following inspection. 
Included is also a penalty for selecting a helicopter 
previously assigned a mission during the model's time frame. 

A reformulation of the short term helicopter - mission 
assignment problem as a network flow model is possible. 
However, due to the success of the model described above (see 
computational results in Chapter IV) and the potential size of 
the network structure, this approach is not further 


investigated. 


E. USER DESCRIPTION, PLANNING SESSIONS AND INTERFACE DESIGN 

To fully understand the role of the two optimization 
mođels, it is necessary to present the context of the 
appropriate planning environment including a description of 
the users. The yearly model is primarily employed once a year 
by a well educated user and therefore can be used as 
demonstrated in chapter IV without an additional interface. 
The applicability of the short term model to assist with 
mission assignment is also demonstrated in chapter IV. 
Unfortunately, the intended user for the short term model 
lacks substantial familiarity with a personal computer. 
Therefore this section proposes a user interface for the short 
term model to assist with complex input and data manipulation 


needs. 
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1. Yearly Planning System 
a. User Description 
The yearly planning system user should have the 
following training, experience and educational requirements: 
@ Complete familiarity with the overall planning process, 


@ Computer skills including input/output problems and 
handling of textfiles, 


€ Basic skills in linear programming and familiarity with 
applicable software packages. 


The use of the system by the maintenance officer in charge 
(usually with a master's degree in air and space technology or 
mechanical engineering) is recommended, but delegation to a 
computer experienced subordinate is possible. Every 
maintenance and repair battalion also has a "Sé-Offizier", who 
is responsible for data processing and the computer 
facilities. 
b. Description of a Planning Session 

After preparing and screening the planning data 
(event information, the planned monthly flight hours, the 
maintenance facility capacities and the initial status of each 
helicopter), a textfile is filled with the required data ina 
Simple coded form. Preferences and prefixed maintenance 
periods can be included. It is helpful to transformers 
achieved optimization model's solution into the more organized 
weekly schedule (now with C inspection recommendations 


included) for further analysis, using a transformation program 
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as described above. Statistics and a graphical overview can 
be prepared for command reviews. The overall time frame for 
a planning session as described depends on the familiarity 
with the system, the solver, the capacity of the personal 
 присес and the difficulty of the plan. Based on 
computational results (see Chapter IV), approximately two 
working days should be sufficient for multiple iterations of 
data input, computer runs, solution review, and manual 
modification of the schedule. 
2. Short Term Planning System 

In addition to the mission assignment optimization 
model as described above, the proposed system (see Figure 3) 
contains graphical interfaces including an update task screen 
for helicopter status changes, a report secticn that includes 
the maintenance status for each helicopter, and a screen for 
generating the updated short term maintenance plan. 

a. Interface Design 

Interface proposals for each screen can be found in 

Appendix E. These interfaces need a not yet developed 
communications manager for screen control and initiating 
computational intermediate steps in the knowledge base, using 
the database information. 

b. User Description 

The recommended user for this part of the planning 


Вүвгетп is one of the two senior NCO’s in the technical command 
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Short Term Planning System 


Status Changes | Ë 


Report 
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Short Term Plan D 


, Planning Data Р i 
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Figure 3 


center. These are usually experienced soldiers with education 
as aircraft technicians and training in maintenance planning 
operations. Their computer skills can be basic. A practical 
briefing of a few days duration should be sufficient to enable 
this group of users to operate the subsystem. Expertise will 
be developed quickly during daily employment. The system 
design should include safety features against wrong input and 
confirmation procedures to prevent unintentional utilization. 
A simple interface structure, easy accessible, controllable 
and with convenient operation possibilities will ensure proper 


use of the system. 
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c. Description of a Planning Session 
(1) Update Process. Every time the technical 
command center gets a vocal (telephone) or written report of 
helicopter status changes, an update screen allows the 
revision for the selected helicopter. These changes include: 
® Failure of an operationally ready helicopter, 
e Category of the failure, 


e Estimated completion times for repairs or maintenance 
measures, 


e Flight hour consumption, 
® Equipment change overs, 
e Special operation limitations, 


ө Renewed availability (mission or maintenance measurement 
completion). 


The helicopter in question, identified by number, will be 
selected and a menu of instances appears. Each instance can 
be selected and changed within a range of offered 
possibilities. The confirmation of the changes concludes the 
update operation. An immediate report generation should be 
included. 

(2) Mission Assignment Process. After determining 
the time frame for the missions in question (from one day to 
a week), a table asks for the mission requirements separated 
by a mission identification number. Simultaneous missions are 
marked with a common group index. Other needed data include 


flight hour requirements, estimated start and landing times, 


px 


description of the mission (combat training program or mission 
order number), spare helicopter requirement, equipment 
requirement and flight restrictions. The short term 
optimization model gives an assignment proposal for each 
mission and spare requirement. After reviewing the result, 
the user can accept the optimization proposal (in part or 
whole), rerun the optimization model, or manually edit the 
selection. The process finishes with a printed result (on 
screen or hardcopy) of the sortie plan. 

(3) Reports and Statistics. A combined status 
report for a particular helicopter or the whole fleet can be 
selected at any time. This feature helps the supervisor in 
his control functions and generates required summaries for the 
next higher command levels. The statistics include: 

€ Percentage of available helicopters, 

€ Percentage of helicopters on mission, 

è Percentage of helicopters in unplanned maintenance, 

è Percentage of helicopters in planned maintenance, 

e Flight hour reserve, 

® Flight hour consumption so far for month and year, 

è Available flight hours for rest of month and year. 
Simple selection of the offered choices should be sufficient 


for generating the report and the statistics. 
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(4) Short Term Maintenance Plan. This special 
screen enables the user to issue his orders to the maintenance 
and repair level facilities. An updated extract of the yearly 
maintenance schedule is generated by selecting helicopters out 
of a candidate list, which includes helicopters already 
scheduled for inspection and those that will soon need 
inspection. The screen allows for comparison between the 


yearly plan and the updated short term version. 
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IV. COMPUTATIONAL EXPERIENCE 


A. THE YEARLY MODEL 

Two data sets (1989 and 1991) taken from the Army Aviation 
Regiment 30 in Niederstetten Air Base (Germany) show the 
applicability of the yearly model. The 1989 data (see Table 
3) requires only two major events and is representative of a 
modestly difficult planning year. The 1991 data (see Table 4) 
pictures a difficult planning year consisting of four 
demanding events and a nonuniform yearly flight hour program. 
Unfortunately, the initial status of each helicopter (hours 
remaining until the next inspection and inspection type) is 
not known for either data set. Realistic estimates are based 
on the fleet’s flight hour reserve, the percentage of 
operationally ready helicopters, and the author’s personal 
experience. All computational results are obtained using an 
ІВМ compatible 486/33 Personal Computer. The model is 
generated using GAMS (Brooke, Kendrick and Meeraus, 1992) and 
solved using ZOOM (Brooke, Kendrick and Meeraus, 1992), XA 
(1993) and OSL (1991). 

The GAMS code, which implements the yearly model can be 
found in Appendix A. The model reports approximately: 

€ 1,200 constraints in 10 blocks, 


Ф 2,600 variables in 18 blocks (600 discrete and 2,000 
continuous variables), 


5ο 


e 9,000 non zero elements. 


Attempts to solve realistic instances of the yearly model 
Optimally result in significant computational effort. 
Heuristic solution procedures were therefore investigated and 


found to provide quality solutions using substantially reduced 


TABLE 9 


DATA SET 1989 


Planned Flight Hours Event 











Helicopter/Hours 
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TABLE 4 


DATA 5ЕТ 1991 











Planned Flight Hours Event 


Helicopter/Hours 





ee 
— 
ee ү >= 


computation time. Two heuristic procedures are used, "LP- 






Rounding" and "Cascade". Both solution procedures rely on the 
relaxed integer solution being a valid information transfer 
tool for the pure integer task. 

B LP-Rounding: After solving the linear programming 
relaxation, binary variables with values near 
zero or one are rounded respectively down or 
up to the nearest integer value and fixed. 


This process can be repeated indefinitely, but 


oe 


Por thrs apmsüe cation two itezations diminish 
the number of decision variables sufficiently 
to allow for a successful solve of the reduced 


mixed integer program. 


B Cascade: A number of possible implementations exist for 
this heuristic technique. For the test problems 
considered, solving twelve (one for each month) 
linear integer programs works best. The first of 
these twelve problems relaxes the integer 
restriction for variables associated with months 
two to twelve. After solving, a new problem is 
generated with month one's variables set equal to 
the values obtained in the previous relaxation, 
month two's variables constrained to be integer, 
and months three to twelve variables relaxed. 
This process repeats until the binary variables 
for all twelve months have been restricted to be 


integer. 


Independent of the heuristic solution procedure, several 
parameter settings are especially important (see Appendix A 
for a complete list). After several test runs uSing various 
parameter values, best results are achieved for: 


BE M  -r Гон om XE,, (flight hourswassigned sto 
helicopter h exceeding regular usage in month t) of 30, 
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® A penadtsy value or ona or as, 

è A penalty value of 0.01 for the elastic variables of 
constraint (6), which provide for the equitable use of the 
helicopter fleet. 

More important than the actual value is the relative magnitude 
of a penalty with respect to the other penalties. 
Specifically, inequable use of the helicopter fleet is viewed 
as an order of magnitude less important than excessive use of 
an “indivitdual@helmeopter “ratio 1/10) 

Solutions from the different solvers (ZOOM, XA,OSL) and the 
two solution methods for the 1991 data set can be compared in 
Appendix B. ZOOM achieves its lowest costs schedules using a 
specified branching order based on the naturel hierarchy of 
the time dependent discrete variables. The solvers XA and OSL 
are applied with the default branch and bound schemes. 

The quality of the solutions does not depend on the solver 
or the heuristic, but ZOOM requires substantially more time 
than XA or OSL, and the Cascade procedure requires slightly 
more time than the LP-Rounding method. All achieved solutions 
exceed the initial relative optimality tolerance values 
(OPTCR) of 0.10 or 0.15 (guaranteeing a solution within 10% or 
15% of optimal respectively). These tolerances apply to each 
individual program run and can therefore not serve as 
computation stopping conditions for the final solution. The 
solutions presented in Appendix B have objective function 
values between 11% and 36% of optimal (using an initial cost 


value of 100). However, the resulting schedules are face- 


34 


valid (i.e. judged implementable by expert opinion) and 
considered superior to manually created plans. 

Table 5 and 6 compare the manual and optimization model's 
yearly schedule for the two data sets. Without any additional 
manual editing, the selected system schedule for 1991 meets 
all requirements (percentage OT operationally ready 
helicopters and flight hour reserve are in the desired ranges) 
and violates the realistic 30 hour planning goal and the 
planned hour per month constraints less than the manual 
solution (fewer planning conflicts). The computed schedule 
for the data set of 1989 shows similar advantages. The 
planned hours per month are now met exactly. Production time 
for both system schedules (including time needed for manual 
adjustment work) constitutes only а small proportion 


(approximately 1/12) of the usual manual completion time. 
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TABLE 5 


COMPARISON SYSTEM SOLUTION - MANUAL ACHIEVED SOLUTIONS 


Abbreviations: Plan - Planned Flight Hours (yearly flight program), 


Man - Manual Planned Flight Hours 


реу! = absolute Deviation from Pl. Hours, Com = Computed Planned Flight Hours 
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TABLE 6 


COMPARISON SYSTEM SOLUTION - MANUAL ACHIEVED SOLUTIONS 


(for Abbreviations see Table 5) 















Simple set 





Solver: XA Solution procedure:  LP-Rounding 
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B. THE SHORT TERM HELICOPTER - MISSION ASSIGNMENT MODEL 

The short term model is tested using data based on past 
experience in mission assignment. A sample test problem 
(including solution) can be found in Tables 7A and 7B. This 
20 mission example represents a normal. one day assignment 
task. Initial input data is listed in Appendix D and XA is 
the solver. The model is also tested successfully with a 
highly complex 40 mission example. This case has demand for 
34 missions with 18 spare requirements and for an unplanned 
exercise requesting six helicopters (five with 45 required 
flight hours each and one with 22 required flight hours). 

The model handles the test cases smoothly when using an 
IBM compatible 486/33 Personal Computer. It produces 
solutions within one to three minutes. All achieved test 
solutions have objective values within five percent of 
optimality. 

The GAMS code implementing the short term model can be 
found in Appendix D. The default branch and bound scheme for 
the XA solver is used. 

A comparison of manual and computed assignments is 
difficult and not appropriate when acknowledging the very 
practical purpose of the system. A program of several months 
in an actual environment should be employed to evaluate the 
utility of the model. Tests with realistic mission sets 
however show consistent and appropriate helicopter selections, 


which pass the "common sense" test. Acknowledged advantages 
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TABLE 7A 


MISSION ASSIGNMENT EXAMPLE 


Mission Required Helo ID Equip. Priority Next Insp. 
Equip. 
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Week of 


Insp 


TABLE 7B 


MISSION ASSIGNMENT EXAMPLE 







- 






Spare helicopter required for mission 


Selected spare helicopter 


of system solutions are the absence of conflicting assignments 
caused by user error (given correct input) and the 


substantially reduced development time. 
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V. CONCLUSION 


This thesis develops two optimization models and proposes 
an interface for a self sufficient PC based decision support 
System as an interactive instrument to construct reliable and 
completely organized heliccpter usage and maintenance plans. 

The yearly maintenance and event planning model produces 
face-valid schedules in substantially less time than current 
manual techniques. The proposed procedure is structured, 
outcomes can be stored for reports, command reviews, and as a 
base for future planning. 

Using the proposed interface to enhance data manipulation 
tasks, the short term planning model for mission assignment is 
able to produce daily to weekly sortie plans in minutes and 
free of user error. 

Using integer linear programming as a foundation for a 
computer based maintenance scheduling and mission assignment 
system can substantially reduce the workload and improve the 


quality of this complex planning process. 
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APPENDIX A IMPLEMENTATION OF YEARLY PLANNING SYSTEM 


A.  LP-ROUNDING METHOD 


STITLE Helicopter Maintenance Scheduling 
SSTITLE A.Sgaslik Thesis 


* Computes a yearly schedule for 45 helicopters with respect 
* to flight hour distribution, D inspections and events. 


* The following program characteristics influence solutions drastically 
and need special attention: 


* 


® The relative termination tolerance OPTCR, which means that GAMS 
will stop and report on the first solution found whose objective 
value is within the specified tolerance of the best possible solution, 

e The scalar DEVWEIGHT, which sets a penalty for no equitable use of the 
helicopter fleet, 

€ The scalar EXPLWEIGHT, which sets a penalty for planning above the 
monthly flight hour guideline for each helicopter, 

e The restriction sec of possible helicopters for event decisions R(T), 

e The upper limits for flight hours for each helicopter and month X.UP 
and XE.UP, 

e The allowed percentage ranges for violating constraints (no penalty 
violations). 


f * + + + * X X X * * + 


| CONTROL OPTIONS S S n ш г 
SOFFUPPER OFFSYMLIST OFFSYMXREF 


OPTIONS 
LIMCOL = LIMROW = 0, SOLPRINT = OFF ,DECIMALS = 2 


0, 
RESLIM = 500000, ITERLIM = 150000, OP'TCR = 0.15, SEED = 3141 
OPTION LP = XA , RMIP - XA , MIP = ХА; 


T month, ) 0512 7 

H helicopter identification / 1*45 / 

E event / 1*4 / 

I iterations for rounding method / 1 / 


ALIAS (T,TP); 


PARAMETERS 
ИЕІСНТІ (Т) penalty for deviations 
WEIGHT2 (T) penalty for deviations 
TOLERA (T) tolerance for no penalty deviation of desired flight hour reserve 
U(T) months without events 
Il 1 
2 1 
3 0 
4 1 
5 0 
6 J 
7 1 
В 1 
9 0 
10 1 
11 B 
12 0 / 
PLHRS(T) planned flight hours for month 
УТ! 526 
2 526 
3 1278 
4 602 
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5 1128 
6 827 

7 526 

8 602 

9 21053 
wO 677 
11 237 
172 743 / 


INHRS (H) initial remaining flight hours for each heiicopter 


/45 0 
44 0 
43 0 
42 2 
41 5 
40 7 
1 10 
2 12 
3 13 
4 15 
5 20 
6 35 
7 45 
8 60 
3 75 
10 80 
ΤΙ 90 
12 95 
13 110 
14 135 
15 145 
16 150 
Τη 155 
18 165 
19 175 
20 185 
39 195 
38 200 
37 205 
36 210 
35 520 
34 225 
33 230 
32 235 
31 240 
30 242 
29 250 
28 260 
27 270 
26 272 
25 278 
24 290 
23 300 
22 300 
21 300 / 


DINSP(H) initial repair level inspection status for each helicopter 
* code description 1-D1 2=D2 3=Dla 4=D3 
/1 


> 
“> 
N = ° Q N P b Ú N F. Q N P &ф WNP PWN S (2 Ñ t-> 
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UJ 
σι 
- P Ú N) к и м ны нә им ыи м ныи 


EVENTHRS(E) event requirements in hours for each helicopter 


yit 100 

2 75 

3 62 

4 85 7 

EVENTNMB (E) event requirements in helicopter numbers 

41 10 

2 8 

3 10 

4 oy 


INSPOK(H,T) one if helicopter - maintenance assignment possible 


EVENTOK(H,E,T) one if helicopter - event - time assignment possible; 


TABLE 
S(T,E) one when event in month 
5-2 3074 
Т 
2 
2,741 
4 
5 J 
6 
7 
8 
Э 1 
10 
is 
12 πω 
TABLE 
R(H,E) one when helicopter event assignment initial possible 
1 2 54 4 
1 1 m 1 
2 1.1545 24 
3 Ὁ, ΥΠ 
4 νυ 
5 1 1! 
6 лт ут 
7 ЗУ 
8 να, τ; 
9 JL 1. 1 
10 о а T 
11 па Д “il 
12 I I I 
13 I D 1 
14 TET 
15 I. 
16 TEL 
T 
18 Ly al: 
19 27 71 1, 
20 Т ШІ 
21 ІТ 1 
22 1 1 
23 ΟΙ 1 
24 a) τ. 
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25 τ. νη 
26 ЖЕРІ” 1 
27 II i 
20 J p 11 
30 Ш 1 1 
BT ТИРЕ І 
32 
33 RENI" n! 
34 NECI I. 
B5 DI 1 
36 EJ] 71 
87 EL Т 
38 met 1 1 
89 H1 1 
40 Ἱ 1 
41 INI 1 1 
42 Er 1 L 
43 pt T 
44 КИЛ «1 
45 / 
SCALARS 
MDAYS average maintenance days for one inspection 
5 7 
DMAXDAYS maximal repair level capacity in days per month 
/ 60 / 
PRODHRS hour production per inspection 
/ 300 / 
OPTSUP desired level of flight hour reserve 


DEVWEIGHT penalty for not equitable fleet usage / 0.01 / 


EXPLWEIGHT penalty for planning above 30 hours limit / 0.1 /; 


* parameter calculations 


OPTSUP = CARD(H) * 150; 
LOOP (T, WEIGHT1 (T) @ КЫ); 
LOOP (T, WEIGHT2 (T) 133 


* Flexible tolerance for deviation from desired flight hour reserve 
ШООР(Т, IF (PLHRS(T) GT (4 * PRODHRS), TOLERA(T) = 3 ); 
IF ( (PLHRS(T) GE (2.5 * PRODHRS)) AND 
(PLHRS(T) LE (4 * PRODHRS)) ,TOLERA(T) 


ЕН; 
IF ( PLHRS(T) LT (2.5 * PRODHRS) , TOLERA(T) 1.1}; 


)3 
INSPOK(H,T) = 1 S((((ORD(T)-1) * 100) GT INHRS(H)) AND 


(ORD(T) GT 1)) ; 

EVENTOK(H,E,T) = 1 $ ( (S(T,E) AND (R(H,E)) AND 

( SUM(TP S ((ORD(TP)) LE (ORD(T))), INSPOK(H,TP)) GT O )) 

OR ( S(T,E) AND (R(H,E)) AND 

( SUM(TP $ ((ORD(TP)) LE (ORD(T))), INSPOK(H,TP)) EQ 0 ) AND 

ШОП) — EVENTHRS(E) - ((((ORD(T)) - 1) “ 100) СТ 0) ) ) ὁ}; 
“------------- MODEL--------------------------------------- 
POSITIVE VARIABLES 

КОНТ) assigned flight hours to helicopter in month 

XE (H, T) assigned flight hours above 30 hour limit 

REMHRS (H, T) flight hours until next inspection 


for helicopter at end of month; 


REMHRS.UP(H,T) = 300; 


* Initialization of remaining hours until next inspection 
LOOP (H, 
REMHRS.FX(H,'0') = INHRS(H)) ; 


* Upper limits for assigned flight hours in month 


1-30; 
ШЕСИРҮН,Т) - 30; 


BINARY VARIABLES 
Z(H,E) one if helicopter is assigned to event e 
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Y(H,T) one if helicopter is assigned maintenance at month t; 


* Fixing preplanned inspections in month 1 
LOOP (H, ΙΕ {ΠΝΗΒΕΒΙΗῚ Εὐ υ, χιπ.” 
ЖСЕКІН 2050 


. 
, 
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POSITIVE VARIABLES 


DEVOPTSP1(T) small penalty deviation from desired flight hour reserve 


DEVOPTSP2(T) deviation from desired flight hour reserve 


DEVNOPEN1(T) allowed deviation from desired flight hour reserve 


ELAST1 (T) elastic variable for equation plan 
ELAST2 (T) elastic variable for equation plan 
ELAST3 (T) elastic variable for equation maxd 


DEVFLEET1(H) elastic variable for equitable fleet usage 


DEVFLEET2(H) elastic variable for equitable fleet usage 


PLNOPEN (T) allowed surplus planning hours in month 

PLSMPEN (T) small penalty eiastic variable for equation plan 
FLNOPEN1 (H) allowed deviation from equitable fleet usaqe 
FLNOPEN2 (H) allowed deviation from equitable fleet usage ; 


* Fixing no and small penalty limits 
DEVOPTSP1.UP (T) 0.25 * OPTSUP; 


DEVNOPEN1.UP(T) = 0.10 * OPTSUP; 

FLNOPEN1.UP(H) = 0.35 * (SUM(T $ (ORD(T) GT 1), PLHRS(T)) / CARD(H)) ; 
FLNOPEN2 .UP (H) - 0.35 * (SUM(T $ (ORD(T) GT 1), PLHRS(T)) / CARD(H)); 
PLNOPEN.UP (T) = 0.05 * PLHRS(T); 

PLSMPEN.UP (T) - 0.10 * PLHRS(T); 


VARIABLE 
COST objective function ; 


EQUATIONS 
OBJ objective function 
PLAN (T) constraint on monthly planned hours 
MAXD (T) repair level capacity constraint 
SMOD (T) lower limit for inspections at month 
DEV1 (T) deviation from desired flight hour reserve constraint 
PET) GH, oi computation of remaining flight hours 
FLSM (H) equitable fleet usage 
NUME (E, T) exact number of helicopters for each event 


ІМ51(Т,Е,Н) no maintenance or assigned flight hours when event 


INS2 (T,H) no assigned hours when maintenance ; 


* minimize 


OBJ.—- COST =E= 100 + SUM(T $ (ORD(T) GT 1), 

WEIGHTI(T) * (PLSMPEN(T) + DEVOPTSP1(T)) ) + 
WEIGHT2(T) * ( DEVOPTSP2(T) + ELAST1(T) 

ELAST2(T) * ELAST3(T)) ) ) 

DEVWEIGHT % SUM(H,DEVFLEET1(H) « DEVFLEET2 (H)) 
EXPLWEIGHT * SUM{ (T,H) 5 {(ORD(T) Ст 17, ЖЕНГЕ. 


+ + + — — 


* subject to 


PLAN (Tigo Окор Ст). 
SUM(H , X(H,T) + XE(H,T) + 
SUM(E $ (S(T,E)) , EVENTHRS(E) * Z(H,E) $ (EVENTOK(H, E, T)))9 
PLHRS(T) - ELAST2(T) + ELAST1(T) + PLNOPEN(T) + PLSMPEN(T) ; 
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=B= 


MAXD(T) Š 


SMOD (T) $ 


DEV1 (T) $ 


ΠΗ) 


REMHRS (H, T) 
SUM(E S 
PRODHRS * 


FLSM(H).. 


NUME (E, T) 


ШОСТ, Е.Н) 


BNSZ (T,H) 


MODEL HELICOPTER 
SOLVE HELICOPTER 


(ORD(T) 
SUM(H, REMHRS (H,T)) 
TOLERA (T) 


ЕЕЕ СТ). 


ΠΡΙ ΟἹ)... 


SISTI) 
(ORD(T) GT 


ΤΕ}, 
ШЕШІНУ) 


SUM(T S 


(SUM( T S$ 


DEVFLEET1 (H) 


Εβδ... 
SUM(H, Z(H, 


ШЕ 5(Т,Е)). 


CT) E 


* iterative solving 


POOP (I, 
БӨОР (Н, 


LOOP (E, 


LOOP (T, 


SOLVE HELICOPTER 


PARAMETERS 


ШЕ 


ТЕ 


SOLVE HELICOPTER 


REPORT1(*, T) 


zE- REMHRS (H,T-1) 


(ORD(T) 
SUM( E ,EVENTHRS(E) * 
(ORD (T) 


ШЕЕ ЕНЕ) СЕЮ. 
ТЕСІГІН, Е) 
ο 5 5(Т,Б), 


ЕБІН, Т) 
LOOP(E 9 S(T,E), 


USING 


SUM(H ,MDAYS * Y(H,T) $ INSPOK(H,T)) 
=L= DMAXDAYS + ELAST3(T) ; 


SUM (H, Y(H,T) $ INSPOK(H,T)) =G= 2; 


=G= OPTSUP - DEVOPTSP1(T) - DEVOPTSP2(T) - 


* DEVNOPEN1(T) ; 


quium. 

- X(H,T) - XE(H,T) - 

EVENTHRS(E) * (Z(H,E) $ (EVENTOK(H,E,T)))) « 
$ INSPOK(H,T)) ; 

GT 1), X(H,T) + XE(H,T) + 

(Z(H,E) S (EVENTOK(H,E,T))))) 
GT 1), PLHRS(T)) / CARD(H) ) « 

- DEVFLEET2(H) « FLNOPEN1(H) - FLNOPEN2(H) ; 


ELE 


E) $ (EVENTOK(H,E, T)) ) =E= EVENTNMB(E) ; 
= ο. ο 5 


(ТІН, Т) 


(ЕУЕМТОК (Н,Е,Т)) + 
$ INSPOK(H,T)) =1= 1; 


ZH, T) 


ΙΤ 30 € SUY TH.T) S INSPOK(H,T)) sLs 1; 


πι; 
USING ΕΜΙΡ MINIMIZING COST; 


SS ӘСЗЕЛІН,Е) 
ο ИЕ ЧЕХ Н ЙЕ} 
ЕН, Е) БО І, 


Т 

0 у; 

Ү.ЕХ(Н,Т) 

X.FX(H, T) 

ΝΠ.) 9 b 

L(H,T) EQ1 , Z.FX(H,E) 
X.FX (H, T) 


ΠΟ. 
ΠΕ ΙΖ 


| nu 
e 


DIO. 
TE τυ. 


USING RMIP MINIMIZING COST; 


MIP MINIMIZING COST; 


= — — — — — — — — — — — — — — — — шы — чш шш ы — ы шш — — — чє =< — — — — 


planned hours for each helo and month; 


зо 


REPORT1 (Н,Т) = X.L(H,T) + XE.L(H,T) + SUM(E $ S(T,E), 
EVENTHRSTE)-* Z.L(H,E)) ; 
REPORT1 (’TOTAL’,T) = SUM (H, X.L(H,T) « XE.L(H,T) + 
SUM{E $ S(T,E), EVENTHRS(E).* Z.L(H,E))); 
REPORT1 ('PLAN ',T) - PLHRS(T); 
PARAMETERS REPORT2(*,T) maintenance decision for helo and month ; 
REBORT2  (Н,Т) S ( ORD(T) GT 1) - Y.L(H,T) ; 
BEBORI2 ('TOTAb',T) $ ( ORD(T) GT 1) = SUM( H, Y.L(H,T)); 
PARAMETERS REPORT3(*,E) event decision for helo and event; 
REPORT3 (H,E) = Z.L(H,E); 
REPORTS (’TOTAL’,E) = SUM(H, Z.L(H,E)); 
PARAMETERS REPORT4(T) summed maintenance day for each month; 
REPORT4 (T) 5 (ORD(T) GT 1) = SUM( H, MDAYS * Y.L(H,T)); 
PARAMETERS REPORT5(*,T) flight hour reserve for each helo and month; 


REPORTS (Н.Т) - КЕМНЕ5.І(Н,Т); 
REPORT5 ('TOTAL',T) - SUM (H, REMHRS.L(H,T)); 
МЕРОКТ5 ('OPTIM',T) - 6750; 
PARAMETERS REPORT6 (T) percentage of ready to fly helicopters; 
REPORT6 (T) $ (ORD(T) GT 1) = 100 - 
(100 * SUM (H $ ((REMHRS.L(H,T) EQ O) OR 
7,1 / CARD(H)); 


EEDIDN-REPORT1:2:1:1; 
DISPLAY REPORT1; 
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OPTION REPORT2:2:1:1; 
DISPLAY REPORT2; 


OPTION REPORT3:2:1:1; 
DISPLAY REPORT3; 


DISPLAY REPORT4; 


OPTION REPORTS:2:1:1; 
DISPLAY REPORTS; 


DISPLAY REPORT6; 


DISPLAY PLSMPEN.L,ELAST1.L,ELAST2.L; 
DISPLAY ELAST3.L; 

DISPLAY DEVOPTSP1.L,DEVOPTSP2.L; 
DISPLAY DEVFLEET1.DLb,DEVFDEET2.D; 
DISPLAY XE.L; 


* Creating input for Pascal transformation program Mosch 
FILE RES /INPUTFIL.PAS/ ; 
PUT RES; 


PUT "INPUTFILE MODEL HELICOPTER" 2 
РТ /; 
PUT /; 
LOOP (H, 
LOOP (T S(ORD(T) GT 1), PUT REPORT1(H,T):5:0 ); 
PUT / 7; 
РОТ /; 
LOOP (H, 
LOOP (T SIORD(T) GT Ш ЮЕ РШ ШҮ Н Т ОЗОТ): 
PUL 77); 
PUT /; 
PUT CARD(E):2:0 /; 
LOOP (H, 
LOOP (E, PUT Z.L(H,E):3:0 ); 
PUT / ); 
PUT /; 
LOOP (H, PUT DINSP(H):3:0 ; 
PUT / ); 
PUT /; 
PUT /; 
LOOP (H, 
LOOP (T, PUT REMHRS.L(H,T):4:0 ); 
PUT / ); 
PUT /; 


B. CASCADE METHOD 


Remark: Only parts different to A. are listed. 


SETS 
Т woda 012 / 


H helicopter identification / 1*45 / 
E event / 1*4 / 
I iterations / 1*12 / 


VARIABLE 
COST objective function ; 


EQUATIONS 
OBJ objective function 
PLANA (T) constraint on monthly planned hours 
PLANB (T) constraint on monthly planned hours 
MAXDA (T) repair level capacity 
MAXDB (T) repair level capacity 
SMODA (T) lower limit for inspections at month 
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SMODB (T) lower limit for inspections at month 


DEV1 (T) deviation of desired flight hour reserve constraint 
FLT1A(H,T) computation of remaining flight hours 

FLT1B(H, T) computation of remaining flight hours 

FLSM (H) equitable fleet usage 

NUMEA (E, T) exact number of helicopters for each event 

NUMEB (E, T) exact number of helicopters for each event 


INSIA(T,E,H) no maintenance and planned hours when event 
INSIB(T,E,H) no maintenance and planned hours when event 


INS2A (T,H) no planned hours when maintenance 
TNS2B (T, H) no planned hours when maintenance ; 


* minımize 


OBJ.. COST =E= 100 + SUM(T S (ORD(T) GT 1), 
( WEIGHT1(T) * (PLSMPEN(T) + DEVOPTSP1(T)) ) + 
( WEIGHT2 (T) * ( DEVOPTSP2(T) + ELAST1(T) 
+ ELAST2(T) + ELAST3(T)) ) ) 
+ DEVWEIGHT * SUM(H,DEVFLEET1(H) + DEVFLEET2 (H) ) 
ΠΕΊΘΗΙ “ 50М( (Т,Н) $ (ORD(T) GT 1), XE(H,T) ) ; 


* subject to 


PLANA(T) $ ((ORD(T) GT 1) AND (ORD(T) LE A)).. 
Serie, A(H,T) + XE(H,T) + 
SUM(E 5 (S(T,E)),EVENTHRS(E) * Z(H,E) $ (EVENTOK(H,E,T)))) =E= 
PLHRS(T) - ELAST2(T) + ELAST1(T) + PLNOPEN(T) + PLSMPEN(T) ; 


PLANB(T) $ (ORD(T) GT A).. 
Seman, X (H,T) + XE(H,T) + 
ЕШИ(Е > (S(T,E)),EVENTHRS(E) * ZC(H,E) S (EVENTOK(H,E,T)))) =E= 
PLHRS(T) - ELAST2(T) + ELAST1(T) + PLNOPEN(T) + PLSMPEN(T) ; 


MAXDA(T) $ ((ORD(T) GT 1) AND (ORD(T) LE A)).. 
SUM(H ,MDAYS * Y(H,T) $ INSPOK(H,T)) 
=L= DMAXDAYS + ELAST3(T) ; 

MAXDB(T) $ (ORD(T) GT A).. 

SUM(H ,MDAYS * YC(H,T) $ INSPOK(H,T) ) 
=L= DMAXDAYS + ELAST3(T) ; 


SMODA(T) $ ((ORD(T) GT 1) AND (ORD(T) LE A)).. 
SUM (H, Y(H,T) $ INSPOK(H,T)) =G= 2; 
SMODB (T) $ (ORD(T) GT A).. 
SUM (H, YC(H,T) $ INSPOK(H, T)) =G= 2; 


ШЕШЕ) 5 (ОКр(Т) СТ 1).. 
SUM(H, REMHRS(H,T)) =G= OPTSUP - DEVOPTSP1(T) - DEVOPTSP2(T) - 
TOLERA(T) * DEVNOPEN1(T) ; 


FLT1A(H,T) $ ((ORD(T) GT 1) AND (ORD(T) LE A)).. 
REMHRS (H,T) =E= REMHRS(H,T-1) - 
ШЕТ) - ХЕ(Н,Т) - 
БИГЕ (5 (Т,Е)) , EVENTHRS(E) * (Z(H,E) $ (EVENTOK(H,E,T)))) + 
БИКӘ 5 (Ү(Н,Т) S INSPOK(H,T)) ; 
ου απ τ 5 (ORD(T) GT А).. 
REMHRS(H,T) =E= REMHRS(H,T-1) - 
ШЕШІ) - XE(H,T) - 
SUM(E $ (S(T,E)) ,EVENTHRS(E) * (ZC(H,E) 9 (EVENTOK(H,E,T)))) + 
PRODHRS * (YC(H,T) $ INSPOK(H,T)) ; 


ЕМІН). SUM(T 5 ((ORD(T) GT 1) AND (ORD(T) LE A)), X(H,T) + XE(H,T) + 
SUM( E ,EVENTHRS(E) * (Z(H,E) $ (EVENTOK(H,E,T))))) + 
SUM(T $ (ORD(T) GT A), X(H,T) + XE(H,T) + 
SUM( E ,EVENTHRS(E) * (ZC(H,E) $ (EVENTOK(H,E,T))))) =E= 
(SUM( T $ (ORD(T) GT 1), PLHRS(T)) / CARD(H) ) + DEVFLEET1 (H) 
- DEVFLEET2(H) « FLNOPEN1(H) - FLNOPEN2(H) ; 


NUMEA(E,T) $ (S(T,E) AND (ORD(T) LE A) 


) 
SUM(H, Z(H,E) $ (EVENTOK(H,E,T)) ) =E= EVENTNMB(E) ; 
BENEBIE.T) $ (S(T,E) AND (ORD(T) GT A)).. 
SUM(H, ZC(H,E) $ (EVENTOK(H,E,T)) ) =E= EVENTNMB(E) ; 


INSIA(T,E,H) $ (S(T,E) AND (ORD(T) LE A)).. 
Z(H,E) $ (EVENTOK(H,E, T)) : X(H,T) / 30 + (Y(H,T) $ INSPOK(H,T)) =L= 1; 
INSIB(T,E,H) $ (S(T,E) AND (ORD(T) GT A)).. 
ZC(H,E) $ (EVENTOK(H,E, T)) « X(H,T) / 30 « (YC(H, T) $ INSPOK(H,T)) 


H 
ңы 
Н 
н 


ENSZAUQT,H) $ (U(T) AND (ORD(T) LE A)).. 
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ХІН,Т) / 30 + (Y(H,T) $ INSPOK(H,T)) =L= 1; 
INS2B(T,H) $ (U(T) AND (ORD(T) GT A)).. 
X(H,T) / 30 « (YCIH, T) $ INSPORNH D1) es, 


MODEL HELICOPTER / ALL/ ; 
* iterative solving 


LOOP (Т, А = ORD(I) + 1 ; 
SOLVE HELICOPTER USING MIP MINIMIZING COST ; 
LOOP Н 
ООР (Т 5. (0Кр (Т) ТЕ Ау х ЕНЕ ЖӨНІНЕ) 2“ 
ХЕЛЕХІН,Т) « ХЕСБИНТЕУ); 
УЕХАЛ) Е ЕЕ: 
REMHRS.FX(H,T) - REMHRS.L(H.T) ; 
LOOP(E $ S(T.E). А.ЕХІН;ЕТ-- ο αν το... 


Н 
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APPENDIX B  TABULATED SOLUTIONS OF THE YEARLY SCHEDULE 


IBM comp. 486/33 


Sol.procedure: LP-Rounding 
Solver: ZOOM 


OPTCR Obj.Val. Approximate time frame: 









ΠΠ Zoe 5 


1 









Total Number of D-level inspections : 35 


Total XE hours : 258 


91. 


IBM comp. 486/33 


Sol.procedure: Cascade 
Solver: ZOOM 


OPTCR Obj .Val. Approximate time frame: 
111.07 


6880 
6954 
6876 


7474 
r = Ts c 





Total Number of D-level inspections : 33 


Total XE hours : 106 


5 2 


IBM comp. 486/33 


Sol.procedure: LP-Rounding 
Solver: XA 
OPTCR Obj.Val. Approximate time frame: 
115.86 


ἀκ; 
Month 





Total Number of D-level inspections : 29 


ШЕСІ ХЕ hours 15. 
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IBM comp. 486/33 


Sol.procedure: Cascade 
ο e XA 


eS ee ee о 
|. e 1 2... νο. θα. во И 
|a |] ο. ы Єз а еш Жа 
ЕГЕН гы Жы επ |. δια a 
_— 5 | z |. 2 «θα. ο πας 
bas te | Л г т. EE 
tr — | 2-4 өле | r7. o- oe АНАН 
Ls. S z iu mm вооа лаи 
te |) i a o_o ἂν ως 
dao 6075 i. 56 
‚ы | 3 | εἴ | πε | e 8 
ла | 2 | 6 |. 20.0 | за 


Total Number of D-level inspections : 30 










м JW |н» 


Total XE hours : 242 
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IBM comp. 486/33 


Sol.procedure: LP-Rounding 
Solver: OSL 


ОРТСЕ Obj.Val. Approximate time frame: 
[133.35 | 15 mimtes ΘΟ. 


100 
2 
Š 
2 






з s Jj 
E 2 а 2. ευ. 1 
НЕ: .- j| s | o | 
E o T < | δε | o | 
ЖЕ ЕД: -- ss | o | 
ЕСТІ С вв 7а а ао | 
Ж ЕЦ: | |o s | o 
mo x 6158 | επ 0 
2 | o2 01 sas | 22 | o Ἢ 
ΜΝ. Γ᾽; ο | ο |. о 


Total Number of D-level inspections : 29 


NO 







N IN [м [м 


Шаға! ХЕ hours ἽΝ 17 
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IBM comp. 486/33 


Sol.procedure: Cascade 
Solver: OSE 


OPTCR Obja ala Approximate time frame: 
124.40 


"| з | во | 50 | το, 
τα t o T өзе.) вә | И 
| 3 | 2 J| usss | вэ» | И 
| ο | 2 |-өт” |. δα. | ΟΠΗ 
| s | 2 j| 5632 | 59.9 |. ΠΝ 
| 6 | wn | sos | o L P ΩΝ 
| 48 | 8 οἱ 2505 1| 2522 2| NN 
| ss |. aa N e. | Ж ЖИН 
сов CUTS о 
| 10 6137 | 80.0 30 
| nu | 2 | өс | ма | K 
__ ө7 | өл | M 


Total Number of D-level inspections : 29 

















N IN TÈ IN ΙΝ e IN IN IN IN IN ро 


Total XE hours : 146 
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APPENDIX C TRANSFORMATION PROGRAM FOR YEARLY 


A.  PASCAL CODE 


Compiler: Borland Pascal Version 6.0 (1990) 
program Mosch (input,output); 
($r*] 
uses Crt; 
( Author : Achim Sgaslik 


Assignment : Thesis "Planning German Army Helicopter 
Maintenance and Mission Assignment" 


Written 2810710793 
Update 
Objective : translate GAMS generated monthly helicopter 
inspection and event planning into a weekly 
schedule 
(------ constant and type definition part ------------------------- ) 


const MAX - 45; 
( number of helicopters ) 


type Idtype - 1..MAX; 
( identification numbers for the helicopters ) 


type Dtype = 1..4; 
( D inspections identification D1,D2,D1a,D3) 


type Inspidtype = 1..16; 
{ complete inspection cycle, both levels } 


type Hrstype = integer; 
{ flight hours } 


type Monthtype 


! 
к 
= 
№ 


суре Віпагсуре = 0..1; 

type Insparraytype - array(Idtype]h of Dtype; 

type Hrsarraytype - array[Idtype] of Hrstype; 

type Montharraytype - array[Monthtype] of integer; 

type Matrixarraytype - array[Idtype,0..12] of Hrstype; 
type Bimatrixtype - array[Idtype,Monthtype] of Binartype; 


( for initial input translation and later use as update structure ) 
EnesSrtructype = record 


Idfield lade ype; 
Reminsp > Hestype; 
Nextinsp : Inspidtype; 


end; 


type Inputstructype = array[Idtype] of Structype; 


{ schedule subtypes } 
type Dschedtype = record 
Devela : Ρεγρε; 
Idnr : Idtype; 
Occup : boolean; 
end; 
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SCHEDULE 


type Cschedtype - record 
Idnr пасуре; 
Occup boolean; 
end; 
type Eventtype - record 
Nr integer; 
Idnr Idtype; 
Occup boolean; 
Inspinevent integer; 
end; 
( Schedule main types ) 
type DScheduletype = record 
Dschedarray array[1..6] of Dschedtype; 
WorkD integer; 
end; 
type DISchedulestructype = array[1..50] of DScheduletype; 
type CScheduletype = record 
Cschedarray array[1..6] of Cschedtype; 
WorkC integer; 
end; 
type CISchedulestructype = array{1..50] of CScheduletype; 
type EvScheduletype = record 
Eventarray array[1..15] of Eventtype; 
end; 
type EventSchedulestructype = array[1..50] of EvScheduletype; 


{ for the Z event decision variable 
type Inevarraytype array[Idtype,1 


variable declaration part 


= — — — — — — = от — = — Ó = — = —< — — — — — — — — — — — — — — — 


from opt. model) 
..10] of Binartype; 


var Dinsp Insparraytype; initial inspection status D - insp ) 
Remhrs Matrixarraytype; Remhrs [Η,Τ] 
Inithrs Hrsarraytype; Inithrs ІН) 
Y Bimatrixtype; Y [H,T] inspection decision ) 
X Matrixarraytype; X [H,T] « XE [H, T] planned hours } 
SumY Montharraytype; от (т, УИ 
DSched DISchedulestructype; 
ClSched CISchedulestructype; 
C2Sched CISchedulestructype; 
EvSched EventSchedulestructype; 
Input Inputstructype; 
Z Inevarraytype; Z [H,E] event decision } 
Start -2..48; Starting week for scheduling D - insp } 
Abwch integer; Deviation variable 


Inputfile,Eventfile,Outputfile 
K,H,T,KO, RM, O, KB, XM, N, D, NO 
Finished,Eventinway,Eventmarker 


procedure Sort ( var SInput 


var SInithrs 





procedure and function declaration part 


Inputstructype 
Hrsarraytype ); 


text; 
integer; 
boolean; 


—— — — moo ooo 


; SRemhrs Matrixarraytype ; 


( Sorting of inputstructure in order of remhrs } 


then begin 


var Büutfterl Buffer? Structype; 
Butfer3a Bufter4 Hrstype; 
dO integer; 
Smallest integer; 
begin 
for J : 1 to (MAX-1) do begin 
Smallest := J; 
for Q := (J + 1) to MAX do begin 
if SInithrs[Q]) < SInithrs [Smallest] 
Smallest := Q; 
end; 
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end; 


if Smallest J then begin 


> 

Bufferl :z SInput [Smallest]; 
Buffer2 :- SInput[J] ; 
Buffer3 := SInithrs [Smallest]; 
Buffer4 .= SInithrs[J]; 
SInput [Smallest] := Buffer2; 
SInput (J) := Bufferl; 
SInithrs[Smallest] :- Buffer4; 
SInithrs[J] :» Buffer3; 

end; 


end; 
end; (Sort) 


procedure Formtranslate ( FoRemhrs : Matrixarraytype 


. 
, 


var FInput : Inputstructype ; FDinsp |: Insparraytype 


var FInithrs : Hrsarraytype ); 
( translates initial data and fills inputstructure ) 
ΕΗ: integer; 


begin 
for H := 1 to MAX do begin 
with FInput [H] do begin 
Idfield := H; 
FInithrs[H] := FoRemhrs[H,0] ; 


case FDinsp[H] of 
1 : begin 

if FoRemhrs [H,0] «-75 then begin 
Reminsp  :- FoRemhrs(H,0]; 
Nextinsp :- 4; 

end; 

if (FoRemhrs[H,0] » 75) and (FoRemhrs [H, 0] 
then begin 


Reminsp = FoRemhrs([H,0] - 75; 
Nextinsp := 3; 

end; 

if (FoRemhrs[H,0] » 150) and (FoRemhrs(H,0] 
then begin 
Reminsp = FoRemhrs[H,0] - 150; 
Nextinsp :- 2; 

end; 


if (FoRemhrs(H,0] » 225) then begin 


Reminsp = FoRemhrs[H,0] - 225; 
Nextinsp :- 1; 

end; 

end; 
2 : begin 

if (FoRemhrs[H,0] «- 75) then begin 
Reminsp  :- FoRemhrs [H,0]; 
Nextinsp :- 8; 

end; 


if (FoRemhrs[H,0] > 75) and (FoRemhrs[H,0] <= 150) 


then begin 


Reminsp - FoRemhrs[H,0] - 75; 
Nextinsp := 7; 
end; 
if (FoRemhrs[H,0] > 150) and (FoRemhrs[(H,0] 
then begin 
Reminsp = FoRemhrs[H,0) - 150; 
Nextinsp := 6; 
end; 


if FoRemhrs[H,0] »- 226 then begin 


Reminsp z FoRemhrs[H,0] - 225; 
Nextinsp :- 5; 
end; 
end; 
3 : begin 
if (FoRemhrs[H,0] «- 75) then begin 
Reminsp := FoRemhrs(H,0]; 
Nextinsp :- 12; 
end; 
if (FoRemhrs[H,0] » 75) and (FoRemhrs[H,0] <= 150) 
then begin 
Reminsp - FoRemhrsí(H,0] - 75; 
Nextinsp :- 11; 
end; 
if (FoRemhrs[H,0] > 150) and (FoRemhrs(H,0] <= 225) 
then begin 
Reminsp := FoRemhrs[H,0] - 150; 


222, 


. 
, 


Nexcinsp :- 10; 
епа; 
if FoRemhrs[H,0] >= 226 then редіп 


Reminsp = FoRemhrs [H,0] - 225; 
Nextinsp :- 9; 
епа; 
епа; 
4 : begin 


if (FoRemhrs[H,0] <= 75) then begin 


Reminsp = FoRemhrs [H, 0]; 
Nextinsp := 16; 
end; 


if (FoRemhrs[(H,0] » 75) and (FoRemhrs([(H,0] <= 150) 
then begin 


Reminsp - FoRemhrs(H,0] - 75; 
Nextinsp :- 15; 
end; 


if (FoRemhrs[H,0] > 150) апа (FoRemhrs[H,0] <= 225) 
then begin 


Reminsp - FoRemhrs(H,0] - 150; 
Nextinsp := 14; 
end; 
if (FoRemhrsíH,0] » 225) then begin 
Reminsp  :- FoRemhrs(H,0] - 225; 
Nextinsp := 13; 
end; 
end; 
end; 
end; 
end; 
end; {Formtranslate} 
function GetInspLength { GDinspfield : Dtype) : integer; 


begin 
case GDinspfield of 
1 : GetInspLength 
2 : GetInspLength 
3 : GetInspLengch 
4 : GetInspLenqth 
end; 
end; (GetInspLength) 


Uu UJ 4 WwW 


procedure Initialize (var IDSched : DISchedulestructype ; 
var IC1Sched,IC2Sched : CISchedulestructype ; 
var IEvSched : EventSchedulestructype); 


var М,І1,12,13,І4 : integer; 


begin 
for M :» 1 to 50 do begin 
wich IDSched(M] do begin 
WorkD := 0; 
for Il := 1 to 6 ао bsgin 
Dschedarray[I1].Occup 
end; 
end; 
with ІС1Ѕсһеа (м) do begin 
WorkC := 0; 
for I2 := 1 to 6 do begin 
CSchedarray[I2].Occup :- False; 
end; 
end; 
with IC2Sched[M] do begin 
WorkC := 0; 
for 13 := 1 to 6 do begin 
CSchedarray [I3].Occup := False; 
end; 
end; 
with IEvSched[M] do begin 
for 14 := 1 to 1S do begin 
Eventarray[I4].Nr :-2 0; 
Eventarray[I4].Occup :- False; 
Eventarray[I4].Inspinevent :-» 0; 
end; 
end; 
end; 
end; (Initialize) 


н 


False; 
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procedure Exfill (var EX Matrixarraytype F 
var EEvSched Eventschedulestructype ; 
var EY Bimatrixtype ; var ESumY Montharraytype; 
var EZ Inevarraytype ; var EDinsp Insparraytype; 
var ERemhrs Matrixarraytype; 
var EInputfile,EEventfile text у; 


( fills all variables with results from optimization model from inputfile ) 


var H,E,T,I,L,K,Buffer,Nrofevents integer; 
Insp Binartype; 
Eventweek, Eventlength integer; 


begin 

readln (EInputfile); 

readln (EInputfile); 

ШЕН :- 1 to 45 do begin 
for T := 1 to 12 do begin 

read (EInputfile,EX[H,T]); 

end; 
readin 

end; 

readln (EInputfile); 


(EInputfile); 


for H := 1 to MAX do begin 
for T := 1 to 12 do begin 
read (EInpucfile,EY[H,T]): 
end; 
readln (EInputfile); 
end; 
readln (EInputfile); 
Drs 1 to 12 do begin 
Buffer := 0; 
for H := 1 to MAX do begin 
if EY[H,T] = 1 then begin 
Buffer := Buffer + 1; 
end; 
end; 
ESumY[T] := 
end; 


Buffer; 


readin (EInputfile,Nrofevents) ; 
for H := 1 to MAX do begin 
for E := 1 to Nrofevents do begin 
read (EInputfile,EZ(H,E]); 
end; 
readin 
end; 


(EInputfile) ; 


readin (EInputfile) ; 


for H := 1 to MAX do begin 
readln (EInputfile,EDinsp[H]); 
end; 
readln (EInputfile); 
for H :- 1 to MAX do begin 
DESEE 0 to 12 do begin 
read (EInputfile,ERemhrs [H,T]); 
end; 
readln (EInputfile); 
end; 
for H := 1 to MAX do begin 
for E :- 1 to Nrofevents do begin 
if EZ(H,E] = 1 then begin 
for I := 1 to E do begin 
readin (EEventfile); 
end; 
read (EEventfile, Eventweek) ; 
read (EEventfile,Eventlength) ; 
reset (EEventfile); 
for L :- 0 to (Eventlength -1) do begin 
K := 1; 
while EEvSched[Eventweek « L].Eventarray[K].Occup do begin 
K ;= K + 1; 
end; 
EEvSched [Eventweek « L].Eventarray[K].Idnr :- H; 
EEvSched[Eventweek « L].Eventarray[K].Occup :- True; 
EEvSched[Eventweek « L].Eventarray[K].Nr :- E; 
end; 
end; 
end; 


St 


end; 
end; (Exfill) 


procedure Printstatistic (var GOutputfile : text ; GRemhrs : Matrixarraytype; 
GY : Bimatrixtype ); 


var Buffer,GT,GH : integer; 


begin 
writeln (GOutputfile); 
writeln (GOutputfile) ; 
writeln (GOutputfile, ‘Percentage of operationally ready Helicopters’), 
for GT :s l1 to 129d6 Bear 
Buffer :- 0; 
for GH :- 1 to MAX do begin 
if (GRemhrs(GH,GT] = 0) or (GY(GH,GT] = 1) then begin 
Buffer :- Buffer + 1; 
end; 
end; 
writeln (GOutputfile, 
(100 - ((Buffer/MAX) * 100) ΤΠ: ΕΠΙ сз 
end; 
writeln (GOutputfile) ; 
end; (Printstatistic) 


procedure Printhours (var TOutputfile : text ; TX : Matrixarraytype ; 
TSumY : Montharraytype !; 


var GT,GH : integer; 


begin 

writeln (TOutputfile) ; 
writeln (TOutputfile,' Planned Hours per month and helicopter’) ; 
write (TOutputfile, ’ у 
for GT := 1 to 12 do begin 

write (TOutputfile,GT:3,' *)-; 
end; 
writeln (TOutputfile); 
writeln (TOutputfile) ; 
for GH := 1 to MAX do begin 

write (TOutputfile,GH:3,’ ο 
for GT := l to 12 do beqin 
write (TOutputfile,TX[GH,GT]:3,* *); 

end; 

writeln (TOutputfile); 
end; 
writeln (TOutputfile); 
writeln (TOutputfile); 
writeln (Toutputfile,' Number of D - inspections per month '); 
tor GT :- 1 to 12 do begin 

writeln (TOutputfile,TSumY[GT]:3,' inspections in month ',GT); 
end; 
writeln (TOutputfile); 
writeln (TOutputfile); 

end; { Printhours } 


function ScheduleC ( SKM : integer ) : integer; 


{ computes start week for C - inspections w.r.t remaining hours 
to the inspections } 


begin 
if SKM < 10 then begin 
ScheduleC := 0; 
end; 
if SKM in (10..20] then begin 
SchedulecC := 1; 
end; 
if SKM in (21..30] then begin 
ScheduleC := 2; 
end; 
if SKM >= 31 then begin 
ScheduleC := 3; 
end; 
end; {Schedulec} 
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procedure Princinpucstruc (PInput : Inputstructype ; var POutputfile : cext; 
PRemhrs : Matrixarraytype ); 


Var H : integer; 


begin 
writeln (POutputfile,' Helicopterdata '); 
writeln (POutputfile,' Nr ID Remhrs  Nextinsp  init.Remhrs'); 
tor H := 1 to MAX do begin 
with PInpuc [H] do begin 
write (POutputfile,H:3,' "js 
write (POutputfile,Idfield:3,' ’); 
write (POutputfile,Reminsp:6,"' Eje 
write (POutputfile,Nextinsp:7,'  '); 
writeln (POutputfile, PRemhrs PInput[H].Idfield,0]:7); 
end; 
end; 
writeln (POutpucfile) ; 
end; {Printinputstruc} 


procedure Printschedule (PEvSched : EventSchedulestructype ; 
PDSched : DISchedulestructype ; 
PClSched,PC2Sched : CISchedulestructype ; 
var POutputfile : text ); 


τα κ 2 Integer; 


begin 
writeln (POutputfile); 
writeln (POutputfile); 
writeln (POutputfile,'**********»* '); 
writeln (POutputfile,'* Schedule * '); 
пле (POUtputfile,’*******«eeeew г), 
writeln (POutpucfile) 


writeln (POutputfile,' Code for Inspections during an event : ”); 


` ` ` чм, 


writeln (POutputfile,' 1 = C1 prior to a C2'); 
writeln (POutputfile,' 2 - C2 DE 
writeln (POutputfile,' 3 - C1 after a C2 ‘); 


writeln (POutputfile, eee eee eee ee eee ee ie ie ой йй йо йй їй ой ой йй Δ πε! ) ; 


for M := 1 to 50 do begin 


writeln (POutputfile,’Week = ”,М); 
writeln (POutputfile,'D-Schedule'); 
Ker 1; 
while PDSched[M].Dschedarray([K].Occup do begin 
write (POutputfile,'ID Ξ ',PDSched(M).Dschedarray[K].Idnr:3,' ”); 
writeln (POutputfile, 
' with D-Inspection : ',PDSched[M]).Dschedarray[K].Dfield,' m 
K := K + 1; 
end; 


writeln (POutputfile); 

writeln (POutputfile,'Cl1-Schedule'); 

К :- 1; 

while PC1Sched( [M] .Cschedarray [K] .Occup do begin 
writeln (POutputfile, 

"ТО = PCISched[M].Cschedarray[K].Idnr:3,' ”); 

K := K + 1; 

end; 


writeln (POutputfile); 
writeln (POutputfile,'C2-Schedule'); 


Ee 1; 
while PC2Sched[M].Cschedarray[K].Occup do begin 
writeln (POutputfile, 
‘ID = ',РС25сһеа (М) .Сзсһедаггау [К] .Ійпг:3, nne 
K := K + 1; 


end; 


writeln (POutputfile); 

writeln (POutputfile,'Event-Schedule'); 

K := 1; 

while PEvSched [M] .Eventarray [K] .Occup do begin 
write (POutputfile, 


'Eventnumber - ',PEvSched[M].Eventarray([K].Nr,' ”); 
write (POutputfile, 

‘ID s ',PEvSched[M].Eventarray(K].Idnr:3,' ”); 
writeln (POutputfile, 

'Inspection - ',PEvSched[M].Eventarray(K].Inspinevent); 
K := K + 1; 


end; 
writeln (POutputfile); 


writeln (POutputfile, f k k kt ft kt ттт Татта К ΙΙ в Р ) ; 
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end; 
end; {Printschedule} 


procedure Eventcheck (var VEvSched : EventSchedulestructype ; 
VT : Monthtype ; var VEventmarker : boolean ; 
Ind : integer ); 

( checks if and which inspection is due during event month ) 


var KMM,LM : integer; 


begin 
VEventmarker := False; 
for KMM := 0 to 3 do begin 
LM: 
if VEvSched[(T-1)*4+1+KMM] .Eventarray[LM] .Nr > 0 then begin 
while VEvSched [(VT-1) *4+1+KMM] .Eventarray [LM] .Occup do begin 
if VEvSched[(VT-1) *4+1+KMM] .Eventarray[LM].Idnr = H then begin 
VEvSched[(VT-1)*4«1«KMM].Eventarray(LM].Inspinevent :- Ind; 
VEventmarker :- True; 
end; 
LM := LM + 1; 
⁄ end; 
end; 
епа‘ 


end; (Eventcheck) 


procedure Graphicschedule ( var IDSched : DiSchedulestructype ; 
var IC1Sched,IC2Sched : CISchedulestructype ; 
var IEvSched : EventSchedulestructype ; 
var IOutputfile : text ); 


ναι η τα τὸ ле дет: 
Fieldwritten : boolean; 


began 
writeln (IOutputfile); 
writeln (IOutputfile,' Graphical Schedule ”); 
writeln (IOutputfile); 
writeln (IOutputfile, ’ D1 2 D2 3 Dla 4 = D3 ’); 
writeln (IOutputfile,"' CI 5 πο Ε Event А Й 
wriceln (IOutputfile,' Week = quarter of a month on horizontal axis'); 
writeln (IOutputfile,' Helicopter ID on vertical axis SIS 
writeln (IOutputfile); 
write (IOutputfile, ’ ШЕ 
for IM := 1 to 16 do begin 
write (IOutputfile,IM:2,' *'); 
end; 
writeln (IOutputfile); 
writeln (IOutputfile); 
for IH :- 1 to MAX do begin 
write πο ρου ο ο πο, ολ. 


1 


* 


"wow 
"í 1 


for IM := 1 to 16 do begin 
Fieldwritten :- False; 
LM = 1; 


if IEvSched [IM] .Eventarray[LM] .Nr > 0 then begin 
while (IEvSched [IM] .Eventarray[LM] .Occup) do begin 
if IEvSched [IM] .Eventarray [LM] .Idnr = IH then begin 
write (IOutputfile,' E '); 


Fieldwritten :- True; 
end; 
LM : LM * 1; 
end; 
end; 
ЫМ := 1; 


while IDSched([IM].Dschedarray [LM] .Occup do beqin 
if IDSched[IM].Dschedarray[LM].1Idnr = IH then begin 


write (IOutputfile,’ '); 
write (IOutputfile,IDSched[IM].Dschedarray(LM].Dfield:1,' '; 
Fieldwritten :- True; 
end; 
LM :- LM + 1; 
end; 
LM :- 1; 


while IC1Sched [IM] .Cschedarray [LM] .Occup do begin 
if IC1Sched [IM] .Cschedarray [LM] .Idnr = IH then begin 
write (IOutputfilej;' ο ο 
Fieldwritten :- True; 
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end; 
LM := LM + 1; 
end; 


ШЕ 1: 


while IC2Sched[IM].Cschedarray[LM].Occup do begin 
if IC2Sched [IM] .Cschedarray (LM) . Idnr 
write (IOutputfile,' 

True; 


Fieldwritten :- 
end; 
LM := LM + 1; 
end; 


S 


гү; 


if not Fieldwritten then begin 


write (IOutputfile, 
end; 
end; 
writeln (IOutputfile); 
end; 
wrrtein (IOutpucfilel; 
writeln (IOutputfile); 
write (IOutputfile, ' г 
БЕРГІМ :- 17 to 32 do begin 
write (IOutputfile,IM:2,' 
end; 
writeln (IOutputfile); 
for IH := 1 to MAX do begin 
write (IOutputfile,IH:2,' 


, 


Дш: 


DU 


ПОЗИ := 17 to 32 do begin 
Fieldwritten := False; 
EM «5 1; 


') 


, 


- IH then begin 


if IEvSched(IM].Eventarray(LM].Nr » 0 then begin 
while (IEvSched(IM].Eventarray[LM].Occup) do begin 


if IEvSched[IM].Eventarray[LM].Idnr = IH then begin 
write (IOutputfile, ' 


Fieldwritten 
end; 
LM := LM + 1; 
end; 
end; 
БМ: = 1; 


Е 


s= True; 


^Y 


while IDSched(IM].Dschedarray[LM].Occup do begin 
if IDSched[IM].Dschedarray[LM].Idnr 
write (IOucputfile, ’ 


write (IOutputfile, IDSched[IM]) .Dschedarray [LM] .Dfield:1,’ 


Fieldwritten := 


LM := LM + 1; 


LM := 1; 


Ages 


True; 


while ICiSched[IM].Cschedarray [LM]. 
if ІС1Ѕсһеа (ІМ) .Сѕсһедаггау [LM] . 
write (IOutpucfile,"' 

ЖЕЛДЕ; 


Fieldwritten :- 


LM := LM + 1; 


LM :- 1; 


* 


τ; 


while IC2Sched[IM].Cschedarray(LM]. 
if IC2Sched[IM].Cschedarray [LM]. 
write (IOutputfile,' 

True; 


Fieldwritten :- 
end; 
LM := LM + 1; 


$ 


„у: 


if not Fieldwritten then begin 


write (IOutputfile, 
end; 
end; 
writeln (IOutputfile); 
end; 


writeln (IOutputfile); 

writeln (IOutputfile); 

write (IOutputfile,"' £r 

for IM := 33 to 48 do begin 
write (IOutpucfile,IM:2, ’ 

end; 

writeln (IOutputfile); 

for IH := 1 со MAX do begin 
write (IOutputfile,IH:2,' 


, 


"S 


e 


p 


. 
’ 


= IH then begin 


Occup do begin 


Derry 


= IH then begin 


Occup do begin 


ike baba 


= IH then begin 
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2s 


for IM :s 33 to 48 do begin 
Fieldwritten :- False; 
LM := 1; 
TE IEvSched [IM] . Eventarray(LM] .Nr > 0 then begin 
while (IEvSched [IM] .Eventarray[LM] .Occup) do begin 
if IEvSched [IM] .Eventarray[{LM) .Idnr = IH then begin 
write “(ТО рд Ел е ENDS 


Fieldwritten :- True; 
end; 
LM := LM + 1; 
end; 
end; 
ILM := i; 


while IDSched[{IM) .Dschedarray [LM] .Occup do begin 
if IDSched [IM]. Dschedarray [LM] . Idnr = IH then begin 
write (IOutputfile,' ‘); 
write (IOutputfile, IDSched [IM]. Dschedarray[LM].Dfield:1,' '); 
Fieldwritten <= [ruc 


LM := LM + 1; 


ІМ :« 1; 
while IC1Sched[IM].Cschedarray[LM].Occup do begin 
if ICI1Sched[IM].Cschedarray[LM].Idnr = IH then begin 
write (IOutputfile ο... 
Fieldwritten :- True; 
end; 
LM : LM + l; 
end; 


τ" ΞΘ 
чћһі1е ІС25сһеа [1м] .С=спеаавтау ЕМІ Occup do begin 
if IC2Sched[IM].Cschedarray[LM].Idnr = IH then begin 
write (τουερυςΕι τς И Т 
Fleldwritten := True; 


ТМГ ге ТУІ: 


if not Fieldwritten then beqin 
write (IOutputfile,' - ') ; 
end; 
end; 
writeln (IOutputfile); 
end; 
end; (Graphicschedule) 


begin 
assign (Eventfile,‘C:\TP\STUDY\Eventfile.pas’) ; 
assign (Inputfile,’C:\TP\STUDY\Inputfile.pas’) ; 
assign (Outputfile,’C:\TP\STUDY\Outputfile.pas’ ) ; 


rewrite (Outputfile); 
reset (Inputfile); 
reset (Eventfile); 


Initialize (DSehed, CiSched, C2Sched, EvSched) ; 

Exfill (X,EvSched,Y,SumY,Z,Dinsp,Remhrs,Inputfile,Eventfile); 
Formtranslate (Remhrs,Input,Dinsp,Inithrs); 

Sort (Input,Remhrs,Inithrs); 


for K :- 1 to MAX do begin 
H Input [K]-Idfuield; 
for T= 1) te 12) do begin 


{ Scheduling D - inspections } 
if Y[H, T] = 1 then begin 
N := GetInspLength (Dinsp (HE) }; 
if T » 1 then begin 
( initial starting point for scheduling D - inspection 
second week in month prior to required completion ) 


Start :- (T-1)*4 - 1; 

Abwch := 0; 

for O := 0 to 1 do begin 
{ check if event prohibits prior scheduling } 
τ ΕΙ; 


if EvSched[Start « O).Eventarray[(L].Nr > 0 then begin 
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while EvSched{Start + 0] Eventarray[L] .Occup do begin 

if EvSched (Start + 0] .Eventarray[(L] .Idnr = H 
then begin 
Abwch := Abwch + 1; 

end; 
СОСЕ Уи; 

епа; 

епа; 
епа; 


{ no advance scheduling if planned hours 
in previous month » 30 

if ((Abwch - 0) or (Abwch - 1)) and (X[H,T-1] » 30) then begin 
Abwch :- 2; 

end; 


( special case first month ) 
else begin 


Start := -1; 
Abwch := 2; 
end; 
Finished := False; 


while not Finished do begin 


( if workload in week too big or if planned hours in 
required completion month, search for next 
available spot 

if (DSched[Abwch « Start].WorkD » 2) or (XÍH,T] » 0) 
then begin 

case (N « Abwch) of 
Bo 45.) begin 
Abwch := Abwch + 1; 
end; 
6 : begin 
1 check if event prohibits inspection period 
reaching into next month 

Eventinway :- False; 

I 1; 

if EvSched(N + Start + Abwch]. 

Eventarray[L].Nr > 0 then begin 
while (EvSched[N + Start + Abwch]. 
Eventarray [L) .Occup) 
and (L <= 5) do begin 
if (Еу5сһеа(М + Start + Abwch]. 
Eventarray[L].Idnr - H) then begin 
Eventinway :- True; 
end; 
L := L + 1; 
end; 
end; 


1f not Eventinway then begin 
Abwch := Abwch + 1; 


end 
else begin 
Finished := True; 
end; 
end; 
7 : begin 
Eventinway := False; 
Less 


if EvSched[N + Start + Abwch]. 
Eventarray[L].Nr » 0 then begin 
while (EvSched[N + Start + Abwch]. 
Eventarray [L] .Occup) 
and (L <= 5) do begin 
if (EvSched[N + Start + Abwch]. 
Eventarray[L].Idnr = Н) then begin 
Eventinway :- True; 
end; 
L := L + 1; 
end; 
end; 


if not Eventinway then begin 
Abwch :s Abwch + 1; 

end; 

Finished :- True; 


end; 
end; 
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end 
else begin 


Finished :- True; 


end; 


end; 


’ 


[ record inspection in schedule and update workload ) 


for NO := 0 to (N-1) do begin 
KO :- 15 
while (DSched[Start + Abwch + NO] .Dschedarray [KO] .Occup) 


and (KO <= 5} do begin 
κο απ κο + 1; 
if KO >= 5S then begin 
writeln (Outputfile, 


Schedule conflict D inspection for helicopter ',H,' ın month ’,T); 


end; 


fou 
s 
18 


епа 
els 


end 
end; 


end; 
end; 
DSched(Start « Abwch « NO].Dschedarray[KO].Occup :- 
DSched [Start + Abwch + NO] .Dschedarray[KO] .Dfield := 
DSched [Start + Abwch + NO] .Dschedarray[KO] .Idnr := H; 
DSched [Start + Abwch + NO] .WorkD := 
h 


DSched {Start + Abwch + NO] .WorkD + 1; 


, 


pdate initial inspection variable for the case of a 


econd inspection during planning period 
Dinsp[H] - 4 then begin 

Dinsp [H] := 1 

e begin 

Dinsp[H] := (Dinsp[H]) + 1; 


, 


True; 
Dinsp [H] ; 


( Scheduling C1,C2 - Inspections using remaining hours in month 


and 
if (Re 
Eve 

ТЕ 


to next inspection ) 


mhrs[H,T] <= 75) and (Remhrs[H,T-1] > 75) then begin 


ntcheck /EvSched,T,Eventmarker,1); 
not Eventmarker then begin 
XM := Remhrs[H,T-1] - 75; 
KM :- ScheduleC (XM); 
while (C1Sched{(T-1) *4+1+KM] .WorkC >= 2) and (KM <= 
do begin 
KM := KM + 1; 
end; 
Κι. 
while (CiSched[(T-1)*4«1«KM].CSchedarray [KB] .Occup) 
and (KB <= 5) do begin 
KB := KB + 1; 
if KB >= 5 then begin 
writeln (Outputfile, 


' Schedule conflict C1 inspection for helicopter ',H,' in month 


end; 
end; 


CiSched[(T-1)*4«1«KM].CSchedarray[KB].Occup :- True; 


ClSched[(T-1)*4«1«KM].CSchedarray[KB].Idnr :- H; 
C1Sched [(T-1) *4+1+KM] .WorkC ;:= 
C1Sched [ (T-1) *4+1+KM] .WorkC + 1; 


end; 


end; 


2) 


E 1S 


if (Remhrs[H,T] «- 150) and (Remhrs[H,T-1] » 150) then begin 
Eventcheck (EvSched,T,Eventmarker,2); 


if 


^ Schedule conflict C2 inspection for helicopter ',H;* 


not Eventmarker then begin 
XM :- Remhrs[H,T-1] - 150; 
KM :- ScheduleC (XM); 
while (C2Sched[(T-1)*441«KM].WorkC »- 2) and (KM «- 
do begin 
KM : KM + 1; 
end; 
йы де г 
while (C2Sched[(T-1) *4+1+KM] .CSchedarray [KB] .Occup) 
and (KB <= 5) do begin 
КВ := КВ + 1; 
if KB >= 5 Chen begin 
writeln (Outputfile, 


end; 
end; 


2) 


in month D); 


C2Sched[(T-1)*4«1«KM].CSchedarray[KB].Occup :- True; 


C2Sched { {T-1) *4+1+KM] .CSchedarray [KB] .Idnr := H; 
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—— 


C2Sched [(T-1) *4+1+KM] .WorkC := 
C2Sched[ (T-1) *4+1+KM) .WorkC + 1; 
end; 
end; 


if (Remhrs[(H,T) <= 225) and (Remhrs[H,T-1] > 225) then begin 
Eventcheck (EvSched,T,Eventmarker, 3); 
if not Eventmarker then begin 
XM := Remhrs[H,T-1) - 225; 
KM := Schedulec (XM); 
while (C1Sched[(T-1) *4+1+KM] .WorkC >= 2) and 
(KM <= 2) do begin 
KM :- KM * 1; 
end; 
а= 1; 
while (C1Sched[(T-1) *4+1+KM) .CSchedarray [KB) .Occup) 
and (KB «- 5) do begin 
KB = КВ + 1; 
if KB >= S then begin 
writeln (Outputfile, 
' Schedule connflict C1 inspection for helicopter ',H,’ in month ',T); 


end; 
end; 
C1Sched[ (T-1) *4+1+KM) .CSchedarray [KB] .Occup := True; 
C1Sched { (T-1) *4+1+KM] .CSchedarray [KB] .Idnr := н; 


C1Sched[ (T-1) *4+1+KM] .WorkC := 
C1Sched[ (T-1) *4+1+KM) .WorkC + 1; 
end; 
end; 
end; 
end; 
Printinputstruc (Input, Outputfile,Remhrs) ; 
Printschedule (EvSched, DSched,C1Sched,C2Sched, Outputfile) ; 
Printstatistic (Outputfile,Remhrs,Y); 
Printhours (Outputfile,X,SumY); 
Graphicschedule (DSched,C1Sched,C2Sched,EvSched,Outputfile); 
close (Inputfile); 
close (Outputfile); 
close (Eventfile); 
end. 


B. RESULTING OUTPUT 


Data Set 1991, solved with solver OSL and solution 


procedure Cascade; listed partly only; 


Helicopterdata 

Nr order of helicopters w.r.t. remaining hours to next D 
inspection; 

ID = helicopter identification; 


Remhrs = remaining hours to next inspection (both levels) ; 
Nextinsp = code numbers 1..16 for next inspection; 
init.Remhrs = remaining hours to next D inspection; 
Nr ΠΡ Remhrs Nextinsp init.Remhrs 

1 43 0 j 0 

2 44 0 T6 0 

3 45 0 + 0 

4 42 2 8 2 

5 41 5 4 5 

6 40 7 16 7 

y 14 10 4 10 
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8 2 15 8 7e» 

9 3 13 122 13 
10 4 15 16 LS 
ШЕЛ) 5 20 4 20 
LZ 6 oS 8 55 
JS Fi 45 12 45 
14 8 60 16 60 
15 9 75 4 Т.Б 
16 10 5 7 80 
17 ШЕЛІ 15 ШЕТ 90 
18 T2 20 15 95 
19 13 55 ES LELO 
20 14 60 7 1-5 
21 ШЕ 70 11 145 
22 16 75 TS TSO 
28 17 5 D 155 
24 18 T5 6 kes 
25 19 2 10 J 75 
26 20 35 14 135 
27 39 45 10 195 
28 38 50 6 200 
29 ο 55 2 205 
30 36 60 14 210 
δα. 35 0 TO 220 
Sz 34 75 6 225 
DS E 5 1 2:90) 
34 32 10 13 235 
35 ση TS 2) 240 
36 30 17 5 242 
Ey 27 25 1 250 
38 28 35 1.» 260 
39 2 45 E 270 
40 26 47 5 272 
41 2Е 53 AL 278 
42 24 65 J 290 
43 23 75 9 300 
44 22 75 5 300 
45 21 75 300 

* * * * xx xx xx xx 
* Schedule * Remark: week 5 and weeks 10 


* * * * kx kk kk k 


Code for Inspections during an event 


1 = C1 praezr cena pas 
ο τος. 
з = СІ агвет а 62 


хххх ХХХ ХХХХХХК ХХХ Ж ЕХХЕХХХХХХХЖХАХХХ АХ ЖК 


70 


48 omitted; 


Week = 1 


D-Schedule 

Πο - 43 with D-Inspection B 
ID = 44 with D-Inspection 4 
ieee 45 with D-Inspection Ji 


Cl-Schedule 
ШЕ 33 


C2-Schedule 


Event-Schedule 


ХхХхххххххХхХхХхХХХХХХХХХХХХХХХХХХХХХХХХЖХХХХХХХХХХХ Ж 


Week = 2 

D-Schedule 

ШЕ 43 with D-Inspection 3 
ШЕ - 44 with D-Inspection : 4 
E= 45 with D-Inspection : 1 
C1-Schedule 

ШЕЕ 11 

Ір = 12 


C2-Schedule 


Event-Schedule 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


Week = 3 

D-Schedule 

W= 43 with D-Inspection 3 
ID = 44 wlth D-Inspection : 4 
NES 45 with D-Inspection : 1 
C1-Schedule 

w= 31 

πο. - 29 


C2-Schedule 
πα - 19 


Event -Schedule 


* k kk ХХХХХХХ ХХЖХХХХХХХХХХХХХХХХХХХХХХХХХХХХХЖ 


Week = 4 

D-Schedule 

Ір = 44 with D-Inspection : 4 
ШИЮ = 40 with D-Inspection : 4 
MES = 5 with D-Inspection : 1 


p 


C1-Schedule 
C2-Schedule 


Event-Schedule 


ck oco oc ooo kx ck ck ok ck oc ok ok oko ke o ke o ke kx ke Mk kx x M kx kx Mk ke kx € Mk ko ko ko ko ko ko ko k k 


Week = 10 

D-Schedule 

μα - 8 with D-Inspection : 4 
C1-Schedule 

ID = τ. 

ID- =. Fi 


C2-Schedule 
ποτ. 33 


Event-Schedule 


Eventnumber = 1 ID 7 O0 Ines ces io n m 
Eventnumber = 2 1D -- те сл 
Eventnumber = 1 Ір = 24 Inspection - 3 
Eventnumber = 1 Тр 5 “Ineeeeeten — 2 
Eventnumber = Т Юр = 027 Тпересетеп - 3 
Eventnumber = 1) Тр 2550 Inspection =“ 3 
Eventnumber - 1 ID- 31 Inspection = 2 
Eventnumber = 1 ID = 34 Inspection = 2 
Eventnumber = 1 Ір = 55 кс = 2 
Eventnumber = 1 ID = 43 Inspection - 3 


*k k xk xk xxx xx xxx x xxx xx xx xxx xxx xxx x xxx xxx xxx xxx xxx 


Percentage of operationally ready Helicopters 


80.0 % at Month T 
84.4 % at Month p 
84.4 5 at Month B 
86.7 *$ at Month 4 
88.9 s at Month S 
80.0 * at Month 6 
82.2 5 at Month 9 
82.2 ç at Month 8 
84.4 % at Month 9 
80.0 % at Month 10 
84.4 % at Month 11 
86.7 % at Month 12 
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Planned Hours per month 
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6 7 8 
0 QS 
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1З 0 0 
0 0 0 
EE 0 
"06 O TRA 
2 O30 
30 50 O 
30 0 0 
0 θη 2ῃ 
0 0 0 
5 ο 30 
0 0 0 
30 ШИ 30 
25 0 0 
ds 0 0 
ος 0 0 
30 0 0 
0 0 0 
0 0 0 
ο 30 0 
ου. ο oO 
о озо 
30 0 0 
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Ой 35 0 
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20 EE O 0 
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0 ο 20 
30 Oso 
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Graphical Schedule 
ШЕГІ 2 - D2, 3 = Dia, 4 = D3, * = C1, $ = C2, E = Event 
Week = week equivalent (48 for a year) on horizontal axis 
Helicopter Identification on vertical axis 

КОО 5 65 7 8 910 11 12 13 14 15 16 
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APPENDIX D IMPLEMENTATION OF SHORT TERM PLANNING SYSTEM 


STITLE MISSION ASSIGNMENT MODEL 
SSTITLE A.SGASLIK Thesis 


* Selection of helicopters for specified missions 


SOFFUPPER OFFSYMLIST OFFSYMXREF 


OPTIONS 
LIMCOL - 0 , LIMROW = 0 , SOLPRINT = OFF , DECIMALS = 2 
RESLIM - 100000 , ITERLIM - 5000, OPTCR - 0.05 , SEED - 3141 
OPTION LP - XA, RMIP - XA, MIP = XA; 

s DEFINITIONS AND DATA------------------------------------------- 

SETS 


M mission identification /1%20/ 
H helicopter identification /1*45/ 
I simultaneous mission groups /1*6/ ; 


ALIAS (H,HA); 


PARAMETERS 
REMINSP (H) remaining hours to next inspection 
Т 0 
2 0 
3 5 
š 5 
5 5 
6 10 
7 10 
8 10 
3 12 
10 15 
a 20 
12 25 
13 25 
14 30 
15 35 
16 35 
17 35 
18 40 
15 45 
20 45 
21 50 
22 50 
23 56 
24 55 
25 60 
26 60 
27 70 
28 75 
20 7S 
30 74 
31 60 
32 60 
33 40 
34 40 
a5 30 
36 30 
37 30 
38 20 
39 10 
40 5 
41 5 
42 5 
43 0 
44 0 
45 0 / 
AVAIL (H) availability of helicopter for mission 

1 0 
2 0 
3 1 
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№ 
ы» 
O O O P. | p. |. P. O P. O p. p pj B. BÍ. H[. B. B. BB. B. B. B. B. p. EG. O [p [B крз сура сун кък сукна 





RATION(H) operation limitation 
description  0-none, 1=CTP only, 2=VFR only, 3= daylight only 
/ 


ох союусу л> N P 


№ 
о 
ын ыроооооммноошооооооооошшооооооморнооооо 


A9 


40 0 
41 0 
42 0 
43 0 
44 0 
45 0 / 


INSPWEEK(H) week of next inspection out of yearly schedule 


Jl 12 
2 Τα 
3 13 
4 13 
5 14 
6 14 
7 15 
8 15 
9 16 
10 17 
11 17 
Та 17 
13 18 
14 19 

p 15 20 
16 21 
T 22 
18 23 
19 24 
20 25 
21 26 
22 26 
23 26 
24 207 
25 28 
26 29 
27 30 
28 31 
29 32 
30 33 
31 33 
32 30 
33 29 
34 29 
35 28 
36 27 
37 26 
38 25 
39 19 
40 15 
41 15 
42 15 
43 15 
44 13 
45 127 


PRIORITY(H) user given priority for helicopter 
* 0 = lowest 3 highest 
EE 


Pp 
ney 
F2? O O O r2 N Ñ Ñ Uu (Q (G (Q I2 P P F2. рә FO NO FO CO NO CO CO CO OO (ON 
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EQUIPMENT (H) equipment code for each helicopter 
«1511 seats 2=5seats 3=400kg 5=ITL 6=ITR 7=Winch 
Dol 


к o со < O, (n b G 
e 


№ 
ww 
QJ NJ C Ui шы м м м UG UG GJ) GJ CO tQ BN FIO fO N HIP OY UT i Ç) Ñ) -σσ δώ ES: O. Q СН 


М. 


NEXTINSP(H) completion length of next inspection in weeks 


L 0.5 
2 T 
3 3 
4 4 
5 5 
6 0:5 
Е ο 5 
8 oms 
9 0.5 
10 JL 
11 1 
12 3 
13 ο 5 
14 ο 5 
15 1 
16 ons 
17 3 
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EN 


LENGTH (M) length of mission in hours 
/ 


о 


= 
= 
лл о «2 Ox I2. Ul Ñ) > P щл л л (ІЛ 1 b° Ñ) N PO 
eo 


N 


REQUEQU (M) required equipment code for mission 
* 1211 seats 2=5seats 3=400kg S=ITL 6-ІТЕ 7=Winch 
/ 


MO CO —J OY Ul iD Ú) @Ñ) HÀ 


= 
о 
NJ Ú) U B Hi z F2. OY OY Ñ P Ñ) Q) (n Qn Ë Ú) N N t 


wN 


SPAREREQU(M) spare helicopter required for mission 
Z 


ою «J OY Ut i юн 
оннынооо 
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MISSIONLIM(M) not acceptable helicopter operation limitation 
ж 1-СТР опіу τις only 3=Daylight only 4=none 
/ 


1 4 
2 4 
3 4 
4 1 
5 Jl 
6 T 
7/ П 
8 2 
3 2 
WOR? 
EI 3 
ША” 3 
23.3 
I4 3 
ШЫ 11 
16 4 
17 4 
18 4 
оа 
20 4 / 


MISGROUP(M) groups missions which are happening simultaneously 


F 
= 
^^ OA (Л Ui i o uS uS UJ UJ GJ GJ GJ NJ) NJ N ο D!) ҥнҥ 


SCALARS 
CONST1 constant for objective function / 3 
CONST2 constant for objective function / 5 


BINARY VARIABLE 
FM(H,M) one if helicopter is selected for mission 


FS (H,M) one if helicopter is selected as spare for mission ; 
LOOP (M S (SPAREREQU(M) EQ 0), FS.FX(H,M) = 0 ); 
POSITIVE VARIABLES 


PENMULT (H) penalty for assigning one helo for more then one mission ; 


VARIABLE 
COST objective function variable; 
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EQUATIONS 


OBJ objective function equacion 

CHO (M) exactly one helicopter for each mission 

CHOS (M) exactly one spare for each mission when required 
HRS (H) length of mission constraint 

MULT (H) multiple missions constraint 


MULTS(H,I) multiple missions constraint for simultaneous missions; 


OBJ.. COST -E- 100 + 
SUM((H,M) 5 (AVAIL(H)), ((CONSTI - PRIORITY(H)) + 
Q.l * XCONSI2 — NEXTINSPOH)J > 
MAX( (EQUIPMENT!H) - REQUEQU(M)) , (REQUEQU(M) - EQUIPMENT(H)) ) + 
0.1 % ІМӘРИБЕКІН) ) ΕΜΠ); 


SUM((H,M) S (AVAIL(H) AND SPAREREQU (M})}, 

( 0.5 CONSTI TPRIOR TTY VEE CON S T EEN EXEN SIENE) 

MAX( (EQUIPMENT'/H) - REQUEQU(M)) , (REQUEQU(M) - EQUIPMENT (H 
0.1 * INSPWEEK(H) ) * FS(H,M) ) + 

SUM(H, 0.5 * PENMULT(H) ); 


} + 
ЕЕ. 


СНОСИТ ЅОМ (Н $ (АУАІ (Н) АМЮ (ОРЕКАТІОМ (Н) NE MISSIONLIM(M))), 
ЕМІН,М) ) =E= 1; 
СНОЅ (М) $ (ЅРАКЕКЕФОО (М) ).. 50М (Н $ (AVAIL(H) AND (OPERATION (H) 
NE MISSIONLIM(M))), FS(H,M) ) =Е= 1; 
HRS(H) $ (AVAIL(H)).. REMINSP(H) =б= 
SUM(M, LENGTH(M)* (FM(H,M) « FS(H,M))); 
MULTS(H,I) $ (AVAIL(H)).. SUM(M S (MISGROUP(M) EQ ORD(I)), FM(H,M) » 
FS(H,M) ) =L= 1 ; 
MULT(H) $ (AVAIL(H)).. SUM(M, FM(H,M) ) =L= 1 + PENMULT(H) ; 


MODEL MISSION /ALL/; 
SOLVE MISSION USING MIP MINIMIZING COST; 


* intermediate non-integer solution : 
If (MISSION.MODELSTAT EQ 9, MISSION.ITERLIM = 10000; 
MISSION.OPTCR = 0.1; 
SOLVE MISSION USING MIP MINIMIZING COST ); 


* infeasible solution or still intermediate non - integer solution 
If ((MISSION.MODELSTAT EQ 9) OR (MISSION.MODELSTAT EQ 4) OR 
(MISSION.MODELSTAT EQ 10), 
LOOP ІН, PRIORITY(H) = 3); 
MISSION.OPTCR = 0.2; 
SOLVE MISSION USING MIP MINIMIZING COST ); 


DISPLAY FM.L; 
DISPLAY ЕБЕ, 
DISPLAY PENMULT.L; 
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APPENDIX E INTERFACE PROPOSALS 


‘Helicopter Planning System 


| Main Menu 


к= Update Screen 

С” Mission Assignment Screen 
m 

M Report Generation Screen 


B Maintenance Planning Screen 





Figure 4 
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Helicopter Planning System 
Update Screen 
Possible Status Changes for Helicopter 34 : 


E Flight Hour Update 

zz Equipment Change 

πας Availability | i 
PR Planning Priority 


ES Operation Limitation 


eee 











— __ 





Figure 5 
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Helicopter Planning System ` 
Flight Hour Update 








| 


` 


Helicopter 34 . 


Current Flight Hours: 635 hrs 45 min; 
Next Inspection C1 in 35 hrs 45 min; 


___ hours __ minutes at__/__/93 


| 





Figure 6 
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Helicopter Planning System 
Equipment Change 


Helicopter 34 Current Equipment Status: 5 Seats ; 


ШЕН 5 Scat 
ШШШ |! Seats - 





ШШ 400 kg Load 79 
END — internal Tank Left | 


E 





ΚΙ Internal Tank Right 





Figure 7 
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Helicopter Planning System 
Availability Code | 





Helicopter 34 Current Availability : Operational 


ШИШ Failure Fuselage 
EEN Failure Engine 







| Estimated Completion Time : _ _ / _ _ / 93 x 





Figure 8 
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Helicopter Planning System > 
Planning Priority 


Helicopter 34 Current Planning Priority: 3 = high; 


ШШШ O = very low a 
πα |! - low | 
ШШ 2 medium " 








ΚΙ EN! | 


Figure 9 
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Helicopter Planning System 
| Operation Limitation 


— PF ——— MM ИС 


Helicopter 34 Current Limitation: 0 = none ; 
CTP only | 

= VFR only 

daylight only 


WI NI = 
| | dI 








Figure 10 
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Helicopter Planning System 
Mission Assignment Screen 


Pagel  11/11/95-12/11/95 


zr дый... 
Mises em а: ТАР 
7” 


ΕΤ e e| ава 
з |2 |2 Πο παπα 
ре | о > Zt ae 

Τι Γι mme end 
WG. s J s - 


Legend: Lim = notallowed operation limitation code 
Еди = Equipment code Prio = Priority code 





Figure 11 
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Schedule Week : 21 - 22 
Maintenanoe Planning 


Yearly Schedule <= 10 hours left Short Term Schedule 


ЕЕ 
21122) [7 12 [C2|22 
12|С221 7| |І41с12 

121 





Figure 12 
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Helicopter Planning System 


! 


Report Generation Screen | 





ШЕ Гір Hour Info Month Ι 


ШЕҢ Flight Hour Info Year | 
ШИ Гг Hour Reserve 
ШЕҢ Percentage planned Maint. 












ШЕ Availability Total ercentage p 





Figure 13 
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